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EXECUTIVE  SUMMARY 

The  overall  objective  of  the  HAVE  TRIM  project  was  to  validate  the  prediction  of 
flight  control  surface  deflections  from  model  equations  of  motion  in  the  unpublished 
Master’s  Thesis  by  Captain  Gary  D.  Miller.  To  accomplish  this,  a  sensitivity  analysis 
was  performed  of  the  model’s  control  surface  deflection  predictions  to  changes  in 
stability  derivatives.  The  Learjet  25  stability  derivatives  were  updated  to  improve  the 
accuracy  of  the  model  equations  of  motion.  Project  HAVE  TRIM  then  compared  the 
actual  aircraft  trimmed  flight  control  surface  deflections  to  the  results  predicted  by  the 
model  equations  of  motion  incorporating  the  updated  stability  derivatives. 

Included  in  this  document  are  the  test  and  evaluation  procedures,  concepts, 
rationale,  results  and  conclusions  of  the  HAVE  TRIM  project.  The  responsible  test 
organization  was  the  412*  Test  Wing.  The  HAVE  TRIM  Test  Team,  members  of  the 
USAF  Test  Pilot  School  Class  03A,  accomplished  all  test  and  evaluation.  This  project 
was  part  of  the  curriculum  for  both  the  USAF  Test  Pilot  School  and  the  Air  Force 
Institute  of  Technology. 

A  technical  review  board  and  safety  review  board  were  accomplished  per  Air 
Force  Flight  Test  Center  procedures.  Also  prior  to  flight  test,  the.  HAVE  TRIM  team 
determined  the  sensitivity  of  model  results  to  variations  in  the  Learjet  25  stability 
derivatives.  From  this  investigation,  the  most  critical  stability  derivatives  were 
determined.  These  critical  stability  derivatives  were  the  primary  focus  of  the  Learjet  25 
aerodynamic  stability  derivatives  updates. 

Project  HAVE  TRIM  included  ten  flights  (18.5  flight  hours)  in  the  Advanced 
Information  Engineering  Services  (formerly  Veridian)  Learjet  25.  The  flights  occurred  at 
the  Air  Force  Flight  Test  Center,  Edwards  AFB,  CA,  during  8-17  October  2003  with 
support  from  Advanced  Information  Engineering  Services. 

A  portion  of  flight  test  was  used  to  update  the  longitudinal  and  lateral-directional 
stability  derivatives  of  the  Leaijet  25  aircraft.  These  higher  fidelity  stability  derivatives 
were  then  incorporated  into  the  model  equations  of  motion  predictions  for  control  surface 
deflections. 

Finally,  the  HAVE  TRIM  team  compared  control  surface  deflection  predictions  to 
flight  test  data  as  a  function  of  eight  variables.  The  effects  of  aircraft  center  of  gravity, 
velocity,  altitude,  asymmetric  thrust  and  sideslip  angle  were  investigated  in  trimmed, 
straight  and  level  flight.  Flight  path  angle  variations  were  then  used  to  determine 
trimmed  climbing/descending  flight  effects.  Trimmed  turn  conditions  were  investigated 
by  varying  bank  angle  or  load  factor.  Finally,  pitch-pointing  capability  was  examined  by 
changing  the  pitch  rotation  rate. 

All  test  objectives  were  met. 
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INTRODUCTION 


GENERAL 

USAF  Test  Pilot  School  (TPS)  and  the  Air  Force  Institute  of  Technology  (AFIT) 
sponsored  the  HAVE  TRIM  project  as  part  of  a  joint  curriculum  in  support  of  a  Master’s 
degree  thesis.  This  project  was  conducted  under  the  authority  of  the  Commandant, 
USAF  TPS.  Additional  guidance  and  technical  requirements  were  provided  by  AFIT. 

This  thesis  was  based  on  six-degree-of-freedom  flight  control  equations  of  motion 
(EOM)  model  developed  to  predict  trimmed  flight  conditions  for  the  Advanced 
Information  Engineering  Services  (formerly  Veridian)  Leaijet  25.  These  equations  are 
referred  to  in  the  remainder  of  this  report  as  “the  model”  or  the  “model  EOM”.  The 
model  was  incorporated  into  a  Matlab®  routine  that  output  the  predicted  required  aircraft 
control  surface  deflections  (6’ s')  to  achieve  a  specified  orientation  and  flight  path  of  the 
vehicle.  Prior  to  flight  testing,  the  model  was  used  to  determine  the  sensitivity  of  the 
predictions  to  uncertainties  in  the  stability  derivatives.  Then,  flight  test  data  were  used  to 
refine  the  Leaijet-25  stability  derivatives  for  incorporation  into  and  refinement  of  the 
model.  Finally,  flight  test  results  were  compared  to  the  model  predictions. 

BACKGROUND 

Automatic  flight  control  systems  were  a  vital  part  of  many  current  aircraft 
designs.  The  control  system  decreased  pilot  workload  and  increased  flight  safety,  and 
was  increasingly  used  to  stabilize  aircraft  that  otherwise  would  not  be  stable.  Also,  with 
the  increased  use  of  Unmanned  Aerial  Vehicles  (UAVs),  the  automatic  flight  control 
system  was  vital  to  achieve  the  required  performance.  Traditionally,  automatic  flight 
control  systems  have  used  linearized  EOM  to  simplify  aircraft  motion  around  an 
equilibrium  point.  Linearization  methods  have  been  studied  extensively  and 
implemented  on  many  current  aircraft.  To  eliminate  limitations  created  by  linearization, 
the  aircraft  EOM  were  solved  directly  for  the  5’s  necessary  to  achieve  a  desired  flight 
path. 


‘Trimmed  flight”  in  the  HAVE  TRIM  project  was  defined  as  zero  net  moments 
on  the  aircraft,  resulting  in  no  angular  acceleration.  Straight  and  level,  unaccelerated 
flight  (SLUF),  the  typical  definition  of  “trim”,  was  a  subset  of  the  HAVE  TRIM 
definition  oif  trimmed  flight.  In  SLUF  maneuvers,  the  sum  of  the  forces  was  assumed  to 
be  zero,  resulting  in  no  accelerations  on  the  aircraft. 

In  the  HAVE  TRIM  project  the  model  EOM  were  solved  for  6’s  at  trimmed  flight 
conditions  and  compared  to  flight  test  results  of  5’s.  The  ability  to  directly  calculate  5’s 
for  a  given  flight  condition  would  reduce  computational  processing  time  that  would  be 
spent  refining  5’s  in  an  iterative  process  to  set  the  specified  flight  condition. 


’  For  clarity,  the  phrase  “control  surface  deflections”  will  be  replaced  by  the  symbol  6  or  6’s  throughout  the  remainder  of  the  HAVE 
TRIM  TIM.  This  designation  either  refers  to  flight  control  surface  deflections  generically  or  to  the  combination  of  the  elevator  (6e), 
aileron  (6a).  and  rudder  (6i)  surface  deflections.  Symbols  6e,  6a,  and  6c  in  the  HAVE  TRIM  TIM  refer  specifically  to  the  identified 
control  surface  only. 


1 


Edwards  Air  Force  Base 

DECEMBER  2003  Air  Force  Flight  Test  Center 

PROGRAM  CHRONOLOGY 

Flight  testing  was  conducted  from  the  Air  Force  Flight  Test  Center  (AFFTC)  at 
Edwards  AFB  California,  in  the  R-2508  complex,  from  8  to  17  October  2003.  Flight 
testing  was  conducted  by  the  HAVE  TRIM  test  team,  USAF  TPS  Class  03A  (The 
CENTURIONS).  All  testing  was  conducted  under  job  order  number  (JON)  M03C1400. 
The  Responsible  Test  Organization  (RTO)  was  the  412*'’  Test  Wing.  Flight  testing 
consisted  of  ten  Learjet-25  sorties  totaling  18.5  flight  hours. 

TEST  ITEM  DESCRIPTION 

The  test  item  was  the  flight  control  system  of  a  Learjet-25  owned  and  operated  by 
Advanced  Information  Engineering  Services  (formerly  Veridian)  and  registered  in  the 
civilian  experimental  category.  It  was  capable  of  in-flight  simulation  of  the  flight  control 
system  and  flight  dynamics  of  other  aircraft  via  the  Variable  Stability  System  (VSS). 
The  left  cockpit  flight  controls  were  production  standard  reversible  Leaijet-25. 

The  right  set  of  flight  controls  was  mechanically  disconnected  from  the  standard 
Learjet-25  flight  control  system  and  connected  to  the  flight  controls  via  a  fly-by-wire, 
hydraulically  assisted,  irreversible  system  via  the  VSS.  Control  feel  on  the  right  side 
controls  was  variable  through  the  VSS  and  hydraulic  actuators.  Control  force  and 
position  plus  sensor  signals  were  input  into  the  VSS  computer  that  computed  the  5’s 
required  to  simulate  the  flight  dynamics  of  the  aircraft.  Disengagement  of  the  fly-by-wire 
system  was  achieved  by  shutting  off  the  hydraulic  power  to  the  flight  control  surface 
actuators  and  opening  dual  bypasses  around  each  actuator. 

For  the  HAVE  TRIM  project,  the  control  feel  had  nominal  characteristics  and  the 
flight  dynamics  were  that  of  the  unaugmented  Learjet-25.  The  VSS  computer  was 
configured  to  accept  “open-loop”  control  position  information  and  Programmed  Test 
Inputs  (PTIs)  sculpted  by  the  HAVE  TRIM  team. 

TEST  OBJECTIVES 

The  overall  test  objective  was  to  validate  the  HAVE  TRIM  model  through 
comparison  of  predicted  5’s  with  flight  test  measured  5’s  at  specified  trimmed  conditions. 
HAVE  TRIM  project  had  three  primary  objectives  in  support  of  the  overall  objective. 
These  primary  objectives  were: 

-  Evaluate  model  sensitivity  to  aerodynamic  stability  derivative  changes. 

-  Determine  the  aerodynamic  stability  derivatives  of  the  test  aircraft. 

-  Compare  model  5’s  predictions  to  flight  test  5’s  results. 

All  the  test  objectives  were  met. 
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TEST  AND  EVALUATION 


GENERAL 

Prior  to  flight  testing,  mathematical  computations  were  performed  to  determine 
the  sensitivity  of  the  model  equations  of  motion  (EOM)  to  variations  in  the  stability 
derivatives.  This  was  accomplished  prior  to  flight  test  to  provide  additional  focus  for 
data  collection  and  analysis. 

The  longitudinal  and  lateral-directional  aerodynamic  stability  derivatives  used  in 
the  model  EOM  were  determined  from  data  collected  during  flight  test.  Aerodynamic 
stability  derivative  determination  was  a  function  of  five  variables:  dynamic  pressure  (q), 
Mach  number  (M),  angle  of  attack  (a),  sideslip  angle  (p)  and  flight  regime.  The 
aerodynamic  stability  derivatives  were  calculated  for  each  trinmied  flight  condition. 

The  heart  of  the  HAVE  TRIM  project  incorporated  the  results  from  the  sensitivity 
analysis  and  the  newly  determined  test  aircraft  stability  derivatives  into  the  model  EOM. 
The  model  EOM  are  Matlab®  routines  developed  in  Captain  Miller’s  draft  Air  Force 
Institute  of  Technology  (AFTT)  Master’s  Thesis.  This  thesis  was  based  on  a  six-degree- 
of-freedom  flight  control  EOM  model  developed  to  predict  trimmed  flight  conditions. 
These  equations  are  referred  to  in  the  remainder  of  this  report  as  “the  model”  or  the 
“model  EOM”. 

The  HAVE  TRIM  project  compared  the  predicted  control  surface  deflection  (5) 
results  from  the  model  EOM  for  specified  trimmed  flight  conditions  to  the  actual  aircraft 
flight  test  5  results.  Eight  variables  were  adjusted  in  trimmed  flight  conditions  for  data 
collection.  These  eight  variables  were  aircraft  CG,  true  airspeed  (V),  pressure  altitude 
(PA  and/or  h),  asymmetric  thrust  (T),  flight  path  angle  (y),  sideslip  angle  (P),  bank  angle 
(([))  and  pitch  rate  (Q). 

Flight  Test  Data  Capture  And  Processing 

Flight  test  data  were  collected  by  the  on-board  Leaijet  25  data  collection  system. 
Aerodynamic  parameters  (a,  P,  0,  0,  W,  h,  V,  T,  y,  M,  nx,  ny,  nz,  P,  Q,  and  R)  and  6’s 
from  flight  test  were  provided  by  Advanced  Information  Engineering  Services  (formerly 
Veridian)  in  Matlab®  files.  Parameter  Estimation  (PEST)  was  used  to  calculate  stability 
derivatives.  For  details  on  working  with  PEST  and  the  HAVE  TRIM  PEST  lessons 
learned,  refer  to  Appendix  K.  All  other  data  processing  and  results  presentation  were 
done  using  Matlab®. 
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Flight  Envelope  Regimes 

The  Learjet-25  flight  envelope  (see  Appendix  A  for  a  detailed  description  of  the 
flight  envelope)  was  divided  into  the  following  three  flight  regimes  or  regions  listed  in 
priority  order  and  characterized  by  different  combinations  of  q  and  M.  These  regions  are 
further  detailed  in  Appendix  B.  Due  to  data  collection  time  constraints  during  HAVE 
TRIM  execution,  data  were  not  collected  in  Regime  3.  This  did  not  prevent  completion 
of  all  test  objectives. 

1)  Medium  M  (0.29  -  0.57)  and  medium  q  (0.76  psi  -  1.53  psi) 

2)  Low  M  (0.21  -  0.38)  and  low  q  (0.31  psi  -  0.76  psi) 

3)  High  M  (>0.44)  and  high  q  (>1.53  psi) 

The  model  EOM  were  developed  for  Regime  1,  the  priority  (center)  region,  where 
all  model  EOM  assumptions  were  valid.  Data  were  collected  outside  Regime  1  to  test  the 
validity  and  influence  of  the  assumptions  on  the  accuracy  of  the  6  predictions,  the 
priority  for  the  three  regimes  was  set  based  on  the  relative  importance  of  model  EOM  per 
regime  to  the  Capt  Miller’s  Thesis. 

For  each  regime  a  center  test  point  was  defined  at  15,000  ft  PA  and  a  specific 
velocity.  Testing  for  the  model  EOM  evaluation  at  Regime  1  was  centered  on  the 
velocity  of  220  knots  indicated  airspeed  (KIAS)  (approximately  0.44  M).  Flight  testing 
in  Regime  2  was  centered  on  150  KIAS  (approximately  0.30  M),  this  corresponds  to 
approximately  1.3  times  stall  airspeed  (1.3  Vstaii^)- 

The  effects  on  the  5’s  of  each  of  the  eight  variables  mentioned  above  were 
examined  separately,  thus  neglecting  the  combined  influence  of  the  variables.  The  results 
from  the  investigation  in  Regime  2  were  used  to  evaluate  the  capability  of  the  model  to 
predict  the  5’s  in  a  flight  condition  close  to  Vstaii  where  the  model  assumptions  begin  to 
lose  accuracy.  The  evaluation  in  Regime  2  determined  how  flexible  the  model  was  in 
predicting  the  5’s  in  a  non-linear  region.  The  results  of  this  analysis  helped  define  the 
applicable  envelope  of  the  model. 

Independent  Variables 

The  effects  of  the  eight  variables  investigated  started  with  CG,  h,  and  V.  Then  the 
effects  of  the  other  variables  were  examined  one  at  time.  Thus,  a  build  up  approach  was 
followed  for  evaluating  the  influence  of  the  eight  variables  on  the  5’s,  going  from  less 
demanding  maneuvers^  to  more  demanding  maneuvers.  Overall,  CG  effects  were 
examined  between  approximately  14%  and  21%"*.  Altitude  sensitivity  was  investigated 
between  8,800  and  20,000  ft  PA.  Velocity  effects  were  determined  between  135  HAS 
(approximately  1.15  Vstaii)  and  250  KIAS  (0.50  M)  at  15000  ft  PA.  The  effect  of 
asymmetric  thrust  was  examined  by  varying  the  thrust  from  each  engine  up  to  the 


2  The  stall  speed  was  computed  based  on  a  gross  weight  of  13,500  lbs.  All  test  points  in  Region  2  were  based  on  this  speed  If  the 
weight  was  different  when  the  test  points  were  flown,  the  stall  speed  was  re-computed  in  flight.  The  test  points  were  not  changed 
and  were  flown  if  within  test  limits. 

3  The  analysis  of  CG,  h,  and  V  influence  needed  only  a  set  of  SLUF  trim  shots. 

^  This  CG  range  corresponded  to  aircraft  gross  weight  between  10,000  and  15,000  lbs. 
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maximum  asymmetric  thrust  or  the  maximum  p  test  limit  (see  Appendix  A  for  details  on 
sideslip  limits).  The  effects  of  P  were  evaluated  from  zero  to  maximum  sideslip  test  limit 
(see  Appendix  A),  left  and  right.  The  longitudinal  pointing  capability  with  respect  to  the 
rate  of  climb/descent  was  determined  by  varying  6  from  15“  nose  high  to  5“  nose  low. 
The  lateral  pointing  capability  was  investigated  through  trimmed  level  turns  by  varying  (j) 
from  5“  to  60“  (p  in  the  turns  was  approximately  zero;  maximum  load  factor  (g)  2.5). 
Finally,  the  pitch  pointing  capability  of  the  aircraft  was  investigated  by  testing  negative 
pitch  rates  up  to  approximately  5.6  degrees  per  second  (°/sec).  This  range  was  limited  by 
the  minimum  load  factor  test  limit  of  +0.2g. 

Model  EOM  Sensitivity  To  Variations  In  Aerodynamic  Stability  Derivatives 

Prior  to  the  flight  test,  the  sensitivity  of  the  model  EOM  to  the  aerodynamic 
stability  derivative  variations  was  determined.  This  determination  was  also  used  in  the 
evaluation  of  the  aerodynamic  stability  derivatives,  since  a  bad  estimation  of  a  specific 
aerodynamic  stability  derivative  would  be  insignificant  if  the  impact  on  the  model  EOM 
was  small.  Also,  the  analysis  conducted  for  this  objective  helped  separate  the  effects  of 
the  derivatives  inaccuracy  from  the  effects  of  the  model  EOM  inaccuracy. 

Aerodynamic  Stability  Derivatives  Determination 

The  aerodynamic  stability  derivatives  were  calculated  to  update  and  refine  the 
model  EOM.  Both  longitudinal  and  lateral-directional  doublet  PTIs  were  performed  at 
each  trim  condition.  Aircraft  response  was  captured  via  the  variable  stability  system 
(VSS)  resident  on  the  Learjet-25.  The  different  combinations  of  M  and  q  in  both  regimes 
was  representative  of  the  aerodynamic  flow  energy,  and  was  used  as  a  variable  for  the 
derivatives  determination.  The  flow  energy  (or  flight  regime;  they  are  synonymous  in 
this  context),  affected  the  derivatives  by  modifying  the  aerodynamic  flow  around  the  tail. 
At  low  M  and  q  (Regime  2),  the  flow  energy  was  expected  to  be  lower  than  in  Regime  1. 
The  maximum  a  the  aircraft  could  reach,  the  control  surfaces  control  power  (both 
longitudinal  and  lateral/directional),  and  the  damping  (both  longitudinal  and 
lateral/directional)  were  affected  by  the  different  flow  energy  and  different  combinations 
of  q  and  M.  In  each  regime  the  effects  of  q,  M,  a  and  P  were  investigated. 

Parallel  6  And  Aerodynamic  Stability  Derivatives  Testing 

At  the  start  of  project  HAVE  TRIM,  the  only  aerodynamic  stability  derivatives 
available  were  referenced  to  SLUF  at  15,000  ft  PA  and  250  KIAS.  Therefore  the 
aerodynamic  stability  derivatives  data  collection  and  the  trimmed  5  investigation  was 
performed  simultaneously.  These  aerodynamic  stability  parameters  were  incorporated 
into  the  model  EOM.  This  updated  model  was  used  to  calculate  6’s.  The  comparison 
between  predicted  and  flight  test  data  5’s  was  used  to  evaluate  the  model  EOM  accuracy 
and  the  aerodynamic  stability  derivatives  accuracy.  In  order  to  determine  the  source  of 
errors  in  both  the  model  EOM  and  the  stability  derivatives,  the  aerodynamic  stability 
derivatives  accuracy  were  also  evaluated  by  using  the  10  Times  Cramer-Rao  Parameter 
(lOCRP)  and  expert  insight. 


5 


DECEMBER  2003 


Edwards  Air  Force  Base 
Air  Force  Flight  Test  Center 


TEST  OBJECTIVES 

Determine  sensitivity  of  model  EOM  results  to  stability  derivative  variation 
Sensitivity  analysis  of  the  longitudinal  and  lateral-directional  stability  derivatives 
was  performed  using  the  original  baseline  stability  derivatives.  The  sensitivity  of  the 
stability  derivatives  is  summarized  in  Table  1  below.  Each  stability  derivative  was  varied 
up  to  ±40%  or  the  variation  that  resulted  in  a  corresponding  ±1.0°  5,  whichever  occurred 
first.  Variations  in  stability  derivatives  ±0  to  10%  that  resulted  in  a  ±1.0°  5  displayed 
“Extreme  Sensitivity”.  Variations  in  stability  derivatives  ±10  to  20%  that  resulted  in  a 
±1.0°  5  displayed  “Significant  Sensitivity”.  Variations  in  stability  derivatives  ±20  to 
30%  that  resulted  in  a  ±1.0°  5  displayed  “Moderate  Sensitivity”.  Variations  in  stability 
derivatives  ±30  to  40%  that  resulted  in  a  ±1.0°  5  displayed  “Slight  Sensitivity”.  Stability 
derivatives  not  listed  in  Table  1  below  did  not  produce  a  ±1.0°  5  with  a  ±40%  variation. 
The  stability  derivatives  in  Table  1  are  listed  in  order  of  decreasing  sensitivity  with  the 
most  sensitive  on  top. 


Extreme  Sensitivity 
(±0  to  10%) 

Significant  Sensitivity 
(±10  to  20%) 

Moderate  Sensitivity 
(±20  to  30%) 

Slight  Sensitivity 
(±30  to  40%) 

Cmo  (SLUF 
and  Climbs) 

Cn5j.  (Sideslip) 

CLa  (SLUF) 

Cma  (SLUF 
and  Climbs) 

Cn5r  (Bank) 

C|p  (Sideslip) 

C|5j.  (Sideslip) 

Cyp  (Bank) 

Ci5a  (Bank) 

C^P  (Bank  and 
Sideslip) 

Cj5a  (Sideslip) 

Cjp  (Bank) 

C|5j.  (Bank) 

Table  1;  Results  of  Stability  Derivatives  Sensitivities 
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Determine  6e  sensitivity  to  stability  derivative  variation. 

Table  2  below  displays  the  longitudinal  stability  derivative  variations  at  a  given 
condition  that  resulted  in  a  5  change  of  ±1°  are  listed  under  the  “Error  (Condition)” 
column  along  with  the  condition  at  which  they  occurred.  Those  stability  derivatives 
whose  variations  of  40%  did  not  result  in  a  5  change  of  ±1°  are  listed  under  the  “Error 
(Condition)”  column  as  “>  40%”. 

SLUF,  sideslip,  climbing,  and  turning  conditions  were  examined  by  varying  the 
stability  derivatives  within  the  Matlab®  routine  ±  40%  and  determining  the  resulting 
change  in  5e.  All  longitudinal  stability  derivatives  described  in  Appendix  J  were 
examined. 


Variable 

Error  (Condition) 

Cmo 

±  10%  (SLUF  and  Climbs) 

>  ±  40%  (Sideslips  and  Turns) 

Clxx 

-  23.8%  (SLUF) 

>  +  40%  (SLUF) 

>  ±  40%  (Climbs,  Sideslips,  and  Turns) 

Cma 

±31%  (SLUF  and  Climbs) 

>  ±  40%  (Sideslips  and  Turns) 

Clo 

>±40%  (All) 

Cdo 

>±40%  (All) 

Cmq 

>±40%  (All) 

Cm8e 

>±40%  (All) 

Table  2:  Sensitivity  of  5e  to  Longitudinal  Stability  Derivatives 


Only  three  stability  derivatives  resulted  in  a  5e  change  that  reached  ±1°  when 
varied  ±40%.  During  SLUF  and  the  climbing  conditions  a  change  in  Cmo  of  only  ±10% 
resulted  in  the  limiting  6e.  This  indicated  an  extreme  sensitivity.  Even  a  slight  error  in 
the  stability  derivate  value  would  produce  a  6e  error  of  greater  than  1°.  During  SLUF,  a 
change  of  Clo  of  -23%  resulted  in  the  limiting  5e.  This  derivative  exhibited  a  moderate 
level  of  sensitivity.  During  SLUF  and  climbing  conditions  a  change  in  Cma  of  ±31% 
resulted  in  the  limiting  8e.  The  derivative  exhibited  a  slight  level  of  sensitivity  with  only 
a  large  error  creating  more  than  1°  of  6e.  CmSe  and  C^q  were  expected  to  be  critical,  but 
were  less  than  slight  sensitivity. 
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Determine  5a  and  5r  sensitivity  to  stability  derivative  variation. 

The  results  are  summarized  in  Table  3  for  the  5a  and  Table  4  for  the  5r.  The 
lateral  directional  stability  derivative  variations  at  a  given  condition  that  resulted  in  a  5 
change  of  ±1°  are  listed  under  the  “Error  (Condition)”  column  along  with  the  condition  at 
which  they  occurred.  Those  stability  derivatives  whose  variations  of  40%  did  not  result 
in  a  5  change  of  ±1°  are  listed  under  the  “Error  (Condition)”  column  as  “>  40%”. 

The  stability  derivatives  listed  in  each  table  have  the  most  sensitive  stability 
derivatives  at  the  top  and  are  listed  in  order  of  decreasing  sensitivity.  SLUF,  sideslip,  and 
turning  conditions  were  examined  for  the  limiting  5a  and  5r.  Stability  derivatives  were 
varied  within  the  Matlab®  routine  ±40%  and  the  resulting  changes  in  5a  and  5^  were 
determined.  All  lateral  directional  stability  derivatives  described  in  Appendix  J  were 
examined. 


Variable 

Error  (Condition) 

C,p 

±  6.9%  (Bank) 

±  17.1%  (Sideslip) 

>±40%  (SLUF) 

Cyp 

-  11.0%  (Bank) 

>  +  40%  (Bank) 

>  ±  40%  (Sideslip  and  SLUF) 

Ci5a 

-  11.1%  (Bank) 

+14.6%  (Bank) 

-  18.7%  (Sideslip) 

+  30.0%  (Sideslip) 

>±40%  (SLUF) 

CnSr 

- 13.2%  (Bank) 

+  18.6%  (Bank) 

-22.1%  (Sideslip) 

>  +  40%  (Sideslip) 

>±40%  (SLUF) 

Cnp 

±  15.4%  (Bank) 

±  28.6%  (Sideslip) 

>±40%  (SLUF) 

Cl5r 

±  15.6%  (Bank) 

±  28.6%  (Sideslip) 

>±40%  (SLUF) 

c,p 

>±40%  (All) 

Cnr 

>±40%  (All) 

Cn5a 

>±40%  (All) 

Cl. 

>±40%  (All) 

C„p 

>±40%  (All) 

Table  3:  Sensitivity  of  5a  to  Lateral-Directional  Derivatives 
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No  stability  derivative  when  varied  ±  40%  produced  a  change  of  1°  of  5a  in 
SLUF.  For  a  banked  condition  Cip  produced  a  6a  change  of  1°  with  only  a  ±  6.9% 
variation.  This  condition  and  stability  derivative  exhibited  extreme  sensitivity.  Even  a 
very  slight  error  in  this  stability  derivative  would  produce  a  5a  error  of  greater  than  1°. 
The  next  five  stability  derivatives  listed  in  Table  3  produced  an  5a  change  of  1°  with  less 
than  a  ±  20%  variation  for  a  banked  condition.  Cip  and  Cipa  produced  a  6a  change  of  1° 
with  less  than  a  ±  20%  variation  for  a  sideslip  condition.  These  stability  derivatives 
indicated  a  significant  level  of  sensitivity,  with  only  slight  errors  in  the  stability 
derivatives  creating  more  than  1°  of  6a  error.  CnSr,  Cnp,  and  Cis,  produced  an  6a  change  of 
1°  with  less  than  a  ±  30%  variation  for  a  sideslip  condition.  These  three  stability 
derivatives  indicated  a  moderate  level  of  sensitivity  during  sideslip.  Moderate  errors  in 
these  stability  derivatives  would  produce  more  than  1°  of  6a.  CnSa  was  expected  to  be 
critical,  but  it  was  not.  Cip,  Cyp,  and  Cnp  were  not  expected  to  be  critical,  but  proved 
otherwise. 
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Variable 

Error  (Condition) 

Cn5r 

±  4.6%  (Bank) 

±  10.4%  (Sideslip) 

>±40%  (SLUF) 

Cyp 

-  4.8%  (Bank) 

+  5.6%  (Bank) 

>  ±  40%  (Sideslip  and  SLUF) 

CnP 

±  5.0%  (Bank) 

±  9.3%  (Sideslip) 

>±40%  (SLUF) 

>±40%  (All) 

Ci5a 

>±40%  (All) 

Ci5r 

>±40%  (All) 

Cip 

>±40%  (All) 

Cnr 

>  ±  40%  (All) 

Cn5a 

>  ±  40%  (All) 

Clr 

>±40%  (All) 

Cnp 

>±40%  (All) 

Table  4:  Sensitivity  of  6r  to  Lateral-Directional  Derivatives 


No  stability  derivative  when  varied  ±40%  produced  a  change  of  1°  of  6r  in  SLUF. 
But,  the  first  three  stability  derivatives  listed  in  Table  4  produced  a  6r  change  of  1°  with 
less  than  a  ±5%  variation  for  a  banked  condition.  Cn&  and  Cnp  also  produced  a  5r  change 
of  1°  with  less  than  a  ±10%  variation  for  a  sideslip  condition.  These  conditions  and 
stability  derivatives  exhibited  an  extreme  sensitivity.  Even  very  slight  errors  in  any  of 
these  stability  derivatives  would  produce  a  6r  error  of  greater  than  1°.  CnSa,  Cipa,  and  Cisr 
were  expected  to  be  critical,  but  they  were  not.  Cyp  and  Cnp  were  not  expected  to  be 
critical,  but  displayed  extreme  sensitivity.  Figures  F34,  F35,  H2,  F38,  F39,  F32,  F33, 
F28,  and  F29  show  the  most  critical  derivatives:  Cip,  Cyp,  Cisa,  CnSr,  Cnp. 

Determine  test  aircraft  stability  derivatives 

Elevator  and  combination  rudder/aileron  doublets  were  performed  after  each  trim 
shot  to  collect  time  history  responses  required  for  computing  the  stability  derivatives. 
The  doublets  were  applied  by  using  a  programmed  test  input  (PTI).  PTI  details  are 
discussed  in  Appendix  M.  This  procedure  was  used  on  all  sorties. 

Appendix  G  summarizes  the  longitudinal  and  lateral-directional  stability 
derivatives  at  each  trim  condition.  All  the  derivatives  were  standardized  to  a  CG  position 
equal  to  22%  of  the  mean  aerodynamic  cord.  Flight  conditions  and  the  corresponding 
derivatives  are  presented.  The  lOCRP  are  also  displayed.  Figures  in  Appendix  F  present, 
for  each  point,  the  95%  confidence  interval,  related  to  the  lOCRP  bounds  calculated  by 
Parameter  Estimation  (PEST).  A  hand  faired  curve  was  drawn  through  the  points  and  the 
equation  of  the  curve  was  presented,  whenever  possible.  The  plots  as  function  of  M,  P 
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and  q  are  presented  separately  for  Regime  1  and  Regime  2.  This  isolated  the  effects  of 
M,  q  and  P  (where  applicable)  from  the  effects  of  the  different  trimmed  5e  in  the  two 
regions.  These  different  trimmed  6e  resulted  from  different  ranges  of  a. 

The  Appendix  G  table  is  derived  from  PEST  output  generator  data  base  page. 
Since  the  settings  for  the  data  base  did  not  permit  display  of  more  than  six  digits  for  each 
number,  some  derivatives  and  lOCRP’s  were  truncated.  This  PEST  limitation  had  the 
greatest  effects  on  stability  derivatives  with  magnitudes  of  0.00010  or  less.  Refer  to 
Appendix  K  -  Working  With  PEST  for  more  details.  Stability  derivatives  with  values 
less  than  0.00010  were  Cpo,  CnSa  and  Ci5r.  In  some  cases  Cisa  and  CnSr  also  had  values  of 
0.00010  or  less. 

Data  were  first  evaluated  by  comparing  the  time  histories  calculated  by  PEST  and 
the  time  histories  collected  onboard.  Figure  FI  shows  an  example  of  a  longitudinal  time 
history  comparison,  and  Figure  F27  shows  an  example  of  a  lateral-directional  time 
history  comparison,  taken  from  the  PEST  output  generator  frame.  The  dotted  lines  were 
calculated  by  PEST  and  the  solid  lines  were  those  collected  from  flight  test.  The 
characteristics  to  examine  when  evaluating  the  quality  of  the  match  between  the  two 
curves  were  the  amplitude  of  the  peaks  after  the  PTI  was  applied  and  the  time  between 
the  peaks.  The  time  histories  matched  very  well.  However,  a  slight  time  delay  was 
observed  in  all  the  a  and  P  time  histories.  PEST  predicted  the  initial  change  in  a  and  p 
after  the  PTI  starting  prior  to  the  initial  change  in  flight  test  data  time  histories.  Since 
PEST  was  able  to  closely  match  all  other  parameters  this  apparent  discrepancy  is  most 
likely  the  result  of  acceleration  effects  on  the  a  and  p  measurment  systems  and/or  the 
incorrect  compensation  for  these  acceration  effects  via  the  onboard  data  acquisition 
system.  Determine  the  source  of  and  correct  for  the  discrepancy  between  PEST 
predictions  of  a  and  p  and  onboard  data  acquisition  values  of  a  and  p  prior  to 
future  stability  derivative  determination  on  the  test  aircraft  (Rl) 

Both  longitudinal  and  lateral-directional  stability  derivatives  were  determined  as  a 
function  of  q  and  M.  Longitudinal  derivatives  were  also  determined  as  a  function  of  a 
and  trimmed  5e.  Lateral-directional  derivatives  were  also  determined  as  a  function  of  p. 
Figures  F2  to  F26  show  respectively  Cma,  Clo,  Clo,  C^q  and  CmSe  as  a  function  of  a,  M 
and  q.  Figures  F28  to  F69  show  respectively  Cnp,  Cnr,  CnSr,  Cip,  Cip  and  CiPa  as  a  function 
of  p,  M  and  q.  The  accuracy  of  the  longitudinal  and  directional  derivatives  were 
satisfactory,  while  the  accuracy  for  the  lateral  derivatives  was  marginal;  especially  Cip 
and  Ci5a.  In  the  following  section  each  plot  will  be  examined. 

Figures  F2,  F3,  and  Figure  HI  show  the  most  critical  longitudinal  stability 
derivatives,  Cma,  Clo  Cmo,  as  a  function  of  a.  The  figures  indicate  that  even  though  these 
derivatives  heavily  affected  the  quality  of  the  predictions  in  the  6e,  they  were  calculated 
with  satisfactory  accuracy  to  minimize  the  effects  of  the  sensitivity  on  the  model  results. 


5  Numerals  preceded  by  any  R  within  parentheses  at  the  end  of  a  paragraph  correspond  to  the  recommendation  numbers  tabulated 
in  the  Conclusions  and  Recommendations  section  of  this  report. 


11 


DECEMBER  2003 


Edwards  Air  Force  Base 
Air  Force  Flight  Test  Center 


Figures  F28,  F29,  F32  thru  F35,  F38  thru  F41,  F44  thru  F47,  F50  thru  F53,  F56 
thru  F59,  F62  thru  F64,  F66,  F68,  F69,  and  HI  show  the  most  critical  lateral  directional 
stability  derivatives,  Cip,  Cyp  Cisa,  CnSr,  Cnp  and  Cjsr.as  functions  of  p,  M,  q,  and  a.  The 
effect  of  these  stability  derivatives  on  the  model  results  are  discussed  below. 

The  profile  drag  coefficient  (Cdo)  plus  the  adverse  yaw  and  the  roll  due  to  rudder 
deflection  coefficients  (CnPa  and  Cipr  respectively)  were  characterized  by  significant 
scatter,  thus  making  their  determination  elusive.  Therefore  plots  of  these  derivatives  are 
not  presented.  As  discussed  in  the  sensitivity  analysis  above,  determination  of  these 
derivatives  was  not  critical  to  the  accuracy  of  the  model.  Determine  the  source  of,  and 
correct  for,  the  signiHcant  data  scatter  of  Cdo,  CnSa  and  Cisr  prior  to  future  attempts 
at  determination  of  these  speciHc  stability  derivatives  (R2). 

Overall,  the  quality  of  the  data  collected,  as  assessed  by  time  history  evaluation 
during  PEST,  was  satisfactory  for  the  longitudinal  and  directional  axis.  The  lateral  data 
quality  was  marginal  using  the  same  evaluation  source. 

Data  quality  was  also  evaluated  by  using  the  lOCRP’s,  which  were  plotted  for 
each  data  point  in  all  figures  of  Appendix  F.  The  lOCRP’s  indicated  the  95%  confidence 
interval  for  each  point.  Overall,  the  lOCRP’s  shown  indicated  most  derivative 
calculation  accuracies  were  satisfactory.  Exceptions  are  further  examined  below. 

As  will  be  detailed  below,  the  quality  of  the  data  was  additionally  evaluated  by 
looking  at  the  5  comparison  results.  In  this  evaluation  the  model  EOM  trimmed  6’s 
calculated  using  the  derivatives  obtained  from  PEST  were  compared  with  the  flight  data. 
The  trimmed  longitudinal  and  lateral  5’s  closely  matched  with  the  model  EOM,  resulting 
in  another  indication  of  satisfactory  accuracy  of  the  longitudinal  and  lateral  derivatives 
calculated.  However,  A  marginal  match  was  found  in  the  lateral  axis.  This  confirmed 
that  the  lateral  derivatives  estimation  presented  some  issues. 

Finally,  the  data  quality  was  assessed  by  soliciting  expert’s  opinions  on  the 
results.  Mr  Chris  Nagy,  NASA  Dryden  aeromodeling  and  PEST  expert,  concurred  with 
the  assessment  of  derivative  calculations  performed  during  project  HAVE  TRIM. 

Therefore,  in  order  to  improve  the  quality  of  the  model  predictions,  either  the 
accuracy  of  the  dihedral  effect  and  aileron  control  power  should  be  improved  or  the 
impact  these  derivatives  on  the  model  should  be  limited. 

Derivatives  as  a  function  of  a 

The  Cma  plot  (Figure  F2)  shows  that  the  pitch  moment  coefficient  was  negative, 
as  was  expected  since  the  airplane  was  speed  stable.  It  also  shows  a  positive  slope, 
indicating  that  as  a  increased  the  aerodynamic  center  of  the  airplane  moved  forward,  thus 
reducing  the  distance  between  the  aerodynamic  center  and  the  CG,  which  reduced  the 
longitudinal  stability  margin.  For  positively  cambered  wings,  like  the  Leaijet-25,  this 
result  matched  the  predictions.  The  lOCRP  indicated  the  results  were  satisfactory. 
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The  Clu  plot  (Figure  F3)  shows  that  the  lift  slope  coefficient  was  positive  and  the 
slope  was  negative.  This  resulted  in  flattening  of  the  lift  curve  as  a  increased  and 
approached  Ostaii.  The  lOCRP’s  indicated  the  data  collected  presented  more  uncertainty 
around  the  center  values  than  the  data  shown  in  the  Cma  plot,  but  were  satisfactory 
because  the  time  histories  matched  well  between  prediction  and  reality. 

The  Clo  plot  (Figure  F4)  shows  that  the  zero  lift  coefficient  slightly  increased 
with  a.  This  characteristic  was  caused  by  lift  produced  by  the  horizontal  tail  that  caused 
the  lift  curve  slope  to  decrease  in  the  region  well  below  stall  where  theory  predicted  a 
constant  slope.  The  lift  curve  was  actually  a  greater  order  curve,  resulting  in  a  slight 
change  of  the  zero  lift  coefficient  with  a.  The  lOCRP  accuracy  of  the  data  collected 
indicated  that  at  high  a  the  derivative  estimation  was  less  accurate  than  at  low  a.  The 
data  were  marginal. 

The  Cmq  plot  (Figure  F5)  shows  that  the  pitch  damping  did  not  change 
significantly  with  a.  The  lOCRP’s  bounds  were  larger  than  for  other  derivatives,  but  Cmq 
was  expected  to  be  hard  to  precisely  determine.  This  data  were  also  considered 
satisfactory. 

The  Cm8e  plot  (Figure  F6)  shows  that  the  elevator  control  power  coefficient  was 
negative,  as  expected.  It  also  shows  that  the  slope  of  the  curve  was  negative.  This 
indicated  as  a  increased,  the  down  wash  produced  by  the  main  wing  incident  on  the  tail 
increased.  This  resulted  in  a  decrease  of  the  local  a  on  the  tail  and  consequently  an 
increase  of  lift  produced  by  the  elevator  (the  elevator  load  is  opposite  the  main  wing  lift). 
This  behavior  matched  theory.  In  addition,  the  accuracy  indicated  by  the  lOCRP’s 
indicated  that  the  data  were  satisfactory. 

Derivatives  as  a  function  of  3 

The  Cnp  and  Cm  plots  (Figures  F28  through  F31)  show  the  weathercock  stability 
and  directional  damping  were  basically  constant  with  the  p.  This  was  expected  since  the 
range  of  P  was  not  large  enough  to  present  a  significant  change  in  the  rudder  lateral  force 
that  could  modify  the  directional  stability.  The  lOCRP’s  indicated  that  the  results  were 
satisfactory. 

The  Cn5r  plots  (Figures  F32  and  F33)  show  that  the  directional  control  power  was 
constant  with  p.  The  quality  of  the  data  was  marginal,  since  the  lOCRP  bounds  showed  a 
significant  level  of  uncertainty. 

The  Cip  plots  (Figures  F34  and  F35)  show  that  the  dihedral  effect  was  essentially 
constant  with  p.  The  quality  of  the  data,  according  to  the  lOCRP’s,  was  satisfactory. 
However,  the  points  in  Regime  2  demonstrate  more  scattering  than  the  points  in  Regime 
1.  It  was  possible  that  the  a  influence  at  low  speeds  affected  the  trimmed  flight  condition 
with  P  more  than  at  higher  speed. 
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The  Cip  and  Ci5a  plots  (Figures  F36  through  F39)  show  that  for  both  Regime  1  and 
Regime  2  the  lateral  damping  and  aileron  control  power  were  nearly  constant  with  p. 

The  quality  of  data  was  satisfactory. 

Derivatives  as  a  function  of  M 

The  Cma,  Clh,  Clo,  Cmq,  and  Cni8e  plots  as  a  function  of  M  are  displayed  in  Figures 
F7  through  F16.  The  Cnp  and  Cnr  plots  as  a  function  of  M  are  displayed  in  Figures  F40 
through  F43-  The  CnSr  and  Cip  plots  as  a  function  of  M  are  displayed  in  Figures  F44 
through  F47.  These  plots  showed  that  for  both  Regime  1  and  Regime  2  the  associated 
stability  derivatives  were  basically  constant  with  the  Mach  number.  This  was  as  expected 
because  the  flight  conditions  were  far  from  the  transonic  region  of  the  envelope. 
Therefore,  there  were  no  effects  on  the  aircraft  stability  derivatives  induced  by  shock 
waves,  the  shift  of  the  aerodynamic  center,  or  the  variation  of  lateral  force  on  the  rudder 
with  M.  The  lift  coefficient  of  an  airfoil  increases  with  M  close  to  the  transonic  region; 
therefore,  the  lateral  force  on  the  rudder  could  be  influenced  by  this  factor.  The  lOCRP’s 
indicated  that  the  results  were  satisfactory. 

The  Cjp  plots  (Figures  F48  and  F49)  show  that  for  Regime  2  that  the  lateral 
damping  was  basically  constant  with  the  M.  However,  for  Regime  1,  the  measurements 
Cip  increased  with  M.  Although  unlikely  since  the  M  range  in  Regime  1  did  not 
approach  a  classically  defined  transonic  region,  the  increase  of  Cip  was  deemed  possible. 

The  Ci6a  plots  (Figures  F50  and  F51)  show  the  same  characteristics  of  Cip,  with  a 
trend  for  the  aileron  control  power,  especially  in  Regime  1,  to  increase  with  M.  The 
theory  does  not  support  the  results;  therefore  further  data  analysis  was  suggested.  One 
possible  explanation  could  be  the  limited  number  of  significant  digits  assigned  by  PEST 
to  the  derivatives.  If  the  numbers  were  not  rounded,  the  plot  would  probably  be  different. 
ReHne  Cisa  as  a  function  of  M  in  Regime  1  during  future  test  aircraft  stability 
derivative  determination  (R3). 

Derivatives  as  a  function  of  q 

The  Cma,  Cta,  Clo,  and  Cmq  plots  as  a  function  of  q  are  displayed  in  Figures  F17 
through  F24.  The  plots  showed  that  for  both  Regime  1  and  Regime  2  the  associated 
stability  derivatives  were  essentially  constant  with  q.  This  was  as  expected  because  the 
flight  conditions  were  far  from  a  high  q  region,  where  structural  loads  may  bend  the 
airplane  and  modify  the  aerodynamic  response.  The  lOCRP’s  indicated  the  results  were 
satisfactory. 

The  Cm8e  plots  (Figures  F26  and  F25)  show  that  for  Regime  2  the  elevator  control 
power  was  nearly  constant  with  q.  However,  for  Regime  1,  the  plot  showed  a  decrease 
with  q,  as  if  the  elevator  was  subject  to  bending  from  the  q.  Although  unlikely,  the 
decrease  of  the  control  power  with  increased  q,  was  possible. 

The  c„p,  Cnr,  Cn8r,  Cip,  Cip  and  CiPa  plots  (Figures  F52  through  F63)  show  that  for 
both  Regime  1  and  Regime  2  that  all  lateral-directional  stability  derivatives  were 
basically  constant  with  the  q.  This  was  expected  since  the  flight  conditions  were  far  from 
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a  high  q  region,  where  the  airplane  might  be  subject  to  structural  deformation.  The  only 
exception  was  in  Regime  2  where  Cib  decreased  with  q.  The  dynamic  pressure  was  not 
expected  to  have  any  influence  on  the  stability  derivatives,  in  Regime  2.  Refine  Cip  as  a 
function  of  q  in  Regime  2  during  future  test  aircraft  stability  derivative 
determination  (R4).  Overall,  the  lOCRP’s  showed  a  relatively  small  95%  confidence 
interval,  indicating  that  the  PEST  estimation  of  the  derivatives  was  reliable. 

Derivatives  as  a  function  of  flight  regime  (indicated  by  a) 

Overall,  the  estimation  accuracy  of  the  longitudinal  derivatives  did  not  present 
significant  differences  between  Regime  1  and  Regime  2.  However,  the  data  were 
characterized  by  more  scatter  in  Regime  2  than  in  Regime  1.  As  highlighted  in  the 
previous  sections,  the  lOCRP’s  were  almost  the  same  between  the  two  regimes,  and  the 
variation  of  the  derivatives  with  a,  M  and  q,  except  for  Cm5e,  was  consistent  and 
predictable.  However,  the  flight  regime  seemed  to  influence  the  absolute  value  of  the 
derivatives  calculated.  If  the  derivatives  were  not  plotted  separately,  an  unexpected 
variation  with  Mach  number  or  dynamic  pressure  would  be  noted,  indicating  that  the 
flight  regime  affected  the  aerodynamic  model.  The  most  important  influence  of  the  flight 
regime  on  the  derivatives  was  that  in  Regime  1  the  a  used  to  evaluate  M  and  q  effects  on 
the  aerodynamic  model  were  significantly  different  than  in  Regime  2.  If  the  maneuvers 
in  Regime  1  and  Regime  2  were  flown  at  same  a,  the  difference  between  the  two  regimes 
would  probably  be  negligible.  The  stabilator  position  was  also  changed  between  Regime 
1  and  Regime  2.  The  position  of  the  stabilator  had  an  effect  on  the  elevator  trim  values, 
and  may  have  changed  the  derivatives  between  the  two  regimes.  Investigate  the  effect 
of  test  aircraft  stabilator  position  on  the  stability  derivatives  (R5). 

Overall,  the  estimated  accuracy  of  the  lateral  derivatives  did  not  present 
significant  differences  between  Regime  1  and  Regime  2.  As  highlighted  in  the  previous 
sections,  the  lOCRP’s  were  almost  the  same  between  the  two  regimes.  The  variation  of 
the  derivatives  with  P,  M  and  q,  except  for  a  few  cases,  was  consistent  and  predictable. 
However,  the  flight  regime  seemed  to  influence  the  absolute  value  of  the  derivatives 
calculated.  If  the  derivatives  were  not  split  into  two  regimes  and  plotted  separately,  an 
unexpected  variation  with  P,  M  or  q  was  noted.  This  indicated  that  unless  the  derivatives 
were  normalized  not  only  for  the  same  CG  position  but  also  for  the  same  flight  regime 
(trimmed  5  and  a  range),  then  the  flight  regime  affected  the  aerodynamic  model.  The 
effect  of  the  stabilator  position  may  have  also  contributed  to  the  discrepancies. 

Figures  F64  through  F69  show  Cnp,  Cnr,  Cip,  Cjp,  CnSr,  and  Ci5a,  each  as  a  function 
of  a.  As  predicted,  the  a  influenced  the  value  of  the  lateral-directional  derivatives.  Cnp 
(Figure  F64)  slightly  increased  with  the  a,  indicating  that  the  blanking  effect  to  the 
rudder  was  compensated  by  a  decrease  in  the  destabilizing  effect  of  the  fuselage.  C„r  and 
Cn8r  (Figures  F65  and  F68)  were  almost  constant  with  the  a,  confirming  that  the  rudder 
was  not  subject  to  significant  blanking  by  the  fuselage  as  a  increased.  The  lateral 
derivatives  seemed  more  heavily  affected  by  the  a;  this  could  explain  the  scatter 
observed  in  the  previous  analysis.  Cip  (Figure  F66)  increased  with  a.  Cip  and  Cisa 
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(Figures  F67  and  F69)  decreased  with  a,  as  expected  since  the  wing  lift  coefficient,  for  a 
given  camber  (6a)  tends  to  decrease  with  a. 

The  a  was  responsible  for  the  different  absolute  value  of  the  derivatives  between 
the  two  regimes.  As  previously  discussed,  a  compensation  for  the  different  a  probably 
would  have  made  the  differences  between  the  two  regimes  negligible. 

Ways  to  isolate  the  M  and  q  effects  from  the  a  effects  include  flying  the  FTT’s  in 
a  wind-up  turn,  so  that  a  could  be  held  in  a  specific  interval  with  minimal  changes  in  M 
or  q.  Another  way  could  be  to  normalize  the  derivatives  not  only  for  CG  position,  but 
also  for  a  and  fly  the  FTT’s  for  evaluating  the  M  and  q  effects  at  the  same  a.  Then, 
normalize  the  longitudinal  derivatives  for  CG  and  for  a  in  order  to  reduce  the  flight 
regime  influence. 

Finally,  the  number  of  lateral  directional  time  histories  usable  for  running  PEST 
was  significaintly  less  than  in  the  longitudinal  case,  mostly  because  the  time  spacing 
between  the  aileron  doublets  and  the  rudder  doublets  was  often  not  short  enough  to 
permit  a  good  estimation  of  the  derivatives^. 

Comparison  of  flight  test  6’s  to  model  predictions  of  8’s  at  trimmed  flight 
conditions 

The  trimmed  6’s  were  determined  at  each  test  point  and  compared  to  the 
predictions  from  the  model  EOM.  The  results  were  plotted  for  each  trim  condition 
outlined  below.  Plots  of  the  results  are  in  Appendix  D.  The  model  6  values  were  plotted 
for  an  average  flight  condition  for  each  test  regime.  The  averaged  flight  condition 
calculated  from  flight  test  points  was  13879  pounds  for  the  aircraft  weight,  12712  feet 
PA,  475  feet  per  second  (ft/s)  true  airspeed  in  Regime  1,  and  315  ft/s  true  airspeed  in 
Regime  2.  Stability  derivatives  used  in  the  model  EOM  were  the  average  values  from 
each  of  Regime  1  and  Regime  2. 

Compare  6(CG)  to  model  predictions. 

Data  from  trim  shots  at  the  center  test  conditions  were  analyzed  to  determine  the 
influence  of  CG  position  on  the  aircraft  aerodynamic  parameters  and  5’s  (see  Test  Matrix 
for  details.  Appendix  C).  The  aircraft  parameters  and  5’s  were  evaluated  in  three 
different  CG  position  subsets.  These  subsets  were  each  characterized  by  a  center  value,  a 
maximum  and  a  minimum.  As  the  aircraft  CG  entered  a  new  subset  during  the  each  test 
sortie,  a  trim  shot  was  performed  at  the  prescribed  flight  condition.  Table  5  defines  the 
CG  subsets. 


CG  subsets 

Subsets 

CG1 

CG2 

CG3 

Max 

21.2 

18.1 

15.3 

Center  value 

19.6 

16.7 

14.7 

Min 

18.1 

15.3 

14 

6  Refer  to  Appendix  K  and  L  for  more  details  regarding  PEST  capabilities  and  using  the  program. 
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Table  5:  CG  Subsets 

Table  6  was  used  by  the  test  team  to  determine  the  CG  position  during  flight,  as  a 
function  of  zero  fuel  weight  (with  four  people  on  board),  total  fuel  weight  and  the  weight 
of  fuel  stored  in  the  fuselage. 


fuselage  weight  (lbs) 

gross  weight  (lbs) 

CG  position 

9570 

5200 

500 

14770 

21.1% 

9570 

1700 

500 

11270 

18.1% 

9570 

1700 

0 

11270 

15.3% 

9570 

1200 

0 

10770 

14.6% 

Table  6:  CG  Determination 


The  CG  position  varied  linearly  with  the  total  fuel  between  5200  lbs  and  1700  lbs. 
Then  the  fuel  stored  in  the  fuselage  was  transferred  to  the  wings.  Due  to  time  and  return 
fuel  constraints,  all  test  points  were  flown  before  the  fuel  was  transferred  to  the  wings. 

The  effect  of  changing  the  aircraft  weight  and  CG  was  compared  to  the 
predictions  from  the  model  EOM.  Since  CG  was  a  linear  function  of  the  aircraft  weight 
for  all  test  points,  the  6’s  were  plotted  as  a  function  of  weight.  Figure  D1  shows  the 
combined  results  for  Regime  1  and  Regime  2. 

For  Regime  1  the  flight  test  results  were  almost  all  within  ±0.5"  of  the  theoretical 
prediction.  However,  the  flight  test  results  did  not  show  a  significant  relationship 
between  5e  and  variations  in  weight;  5e  was  constant.  According  to  theory,  a  3000  pound 
increase  in  weight  would  produce  a  1.0°  decrease  in  5e  in  Regime  1,  but  this  was  not 
observed  during  flight  test. 

In  Regime  2  the  data  were  more  scattered  and  also  did  not  match  the  predicted 
trend.  The  Regime  2  flight  test  results  were  on  average  about  1.0°  higher  than  the 
predicted  value.  Again  the  flight  test  6e  were  nearly  constant  for  all  weights  in  the  test 
envelope.  Theoretically,  a  1500  pound  increase  in  weight  would  produce  a  1.0°  decrease 
in  6e  in  Regime  2.  Overall,  the  impact  of  aircraft  weight  did  not  have  a  measurable  effect 
on  the  5e  value  from  flight  test  data.  Theoretically  the  impact  was  expected  to  be  small, 
but  not  as  small  as  the  actual  data  indicated.  The  significant  difference  between  Regime 
1  and  Regime  2  6e  were  the  result  of  a  change  in  trimmed  stabilator  position. 

The  5a  and  5r  were  predicted  to  be  zero  in  trimmed  flight.  Small  values  of  5a  and 
5r  were  measured  in  flight.  These  deviations  result  from  the  impossibility  to  be  perfectly 
trimmed  with  no  lateral-directional  inputs  during  actual  flight  test.  There  were  no 
significant  values  of  5a  and  5r  in  the  weight  investigation  data  points. 

The  plot  of  a  verses  weight  is  shown  in  Figure  D2.  Both  Regime  1  and  Regime  2 
are  shown  on  the  plot.  Figure  D3  shows  the  effect  of  weight  on  0.  For  Regime  1,  the 
flight  test  data  are  scattered  within  2°  of  the  predicted  curve,  with  all  the  0  data  points 
within  1°  of  the  expected  values.  The  a  values  were  directly  affected  by  the  airspeed  at 
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which  the  test  point  was  flown.  A  10  ft/s  change  in  true  airspeed  resulted  in  an  average 
0.5°  to  1.5°  change  in  a.  The  a  and  6  values  were  more  constant  than  the  predicted 
curve.  The  accuracy  of  the  in-flight  a  and  0  measurements  contributed  to  the  scatter  in 
the  data.  The  a  measurement  varied  ±1.0°  to  2°,  and  the  0  measurement  varied  up  to 
±1.0°.  Therefore,  the  weight  effect  on  these  flight  parameters  was  not  as  large  as 
expected. 

In  Regime  1  the  flight  test  data  were  consistently  1°  to  2°  lower  than  the  predicted 
values.  The  flight  test  data  showed  the  same  scatter  in  Regime  2  as  in  Regime  1,  with  0 
nearly  linear  around  the  average  value  and  a  scattered  more  about  the  average  value. 
The  larger  discrepancy  for  a  and  0  in  Regime  2  may  be  due  to  the  fact  that  the  equations 
were  developed  for  Rl,  and  small  aerodynamic  non-linearities  and  the  changed  position 
of  the  stabilator  could  cause  an  error  in  the  predictions. 

The  test  point  set-up  had  a  major  impact  on  variations  between  the  flight  test  data 
and  the  theoretical  predictions.  The  flight  test  velocity  had  the  largest  impact  on  changes 
in  the  5e  values.  Small  variations  in  test  point  velocity  produced  a  large  change  in  8e 
values.  The  effect  of  W  variations  was  less  significant,  and  the  effect  of  changing  PA 
was  very  small.  The  test  set-up  variations  in  P  and  flight  path  angle  from  the  specified 
test  conditions  were  generally  less  than  ±1.0°  for  all  test  points  and  did  not  significantly 
impact  the  results.  The  bank  angle  varied  up  to  ±8.0°  from  the  specified  test  conditions. 
Often  there  was  a  small  P  (±0.5°/sec),  that  was  not  observed  in  flight,  but  impacted  the 
flight  test  data.  The  (t>  measurement  also  had  the  largest  oscillations  about  the  steady- 
state  value,  up  to  a  ±1.0°.  The  combined  errors  in  (|)  values  had  a  significant  impact  on 
the  lateral-directional  results.  Combining  variations  resulted  in  a  combined  change  in  the 
5  value,  which  accounts  for  some  of  the  scatter  seen  in  the  flight  test  data. 

Overall,  the  effect  of  a  produced  a  larger  change  in  the  measured  6e  than  the 
change  in  aircraft  weight  (and  corresponding  CG).  The  weight  effect  was  masked  by  the 
changes  in  a  from  the  predicted  values.  Airspeed  also  had  a  large  effect  on  6e 
measurements  in  flight.  The  scatter  in  the  data,  especially  in  Regime  2,  could  be 
explained  by  the  variances  in  the  a  and  velocity  at  each  test  point  from  the  ideal 
condition  used  for  the  model  EOM  calculation.  When  the  identical  flight  condition 
parameters  were  used  in  the  model  EOM  the  5  values  closely  matched  the  flight  test  data 
(see  Appendix  E  for  examples). 

Compare  8(V)  to  model  predictions. 

The  effect  of  changing  the  aircraft  velocity  in  flight  test  was  compared  to  the 
theoretical  prediction.  The  results  for  both  Regime  1  and  Regime  2  are  plotted  in  Figure 
D4.  The  flight  test  data  followed  the  same  trend  predicted  by  the  model  EOM,  with  a 
small  bias.  The  data  bias  was  more  pronounced  in  Regime  2. 

In  Regime  1,  the  flight  test  5e  was  consistently  higher  than  the  theoretical 
prediction  by  approximately  0.25°.  All  flight  test  data  were  within  0.0°  to  0.5°  of  the 
model  EOM  prediction.  The  effect  of  velocity  was  not  affected  by  small  changes  in  the 
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aircraft  weight  and  test  point  altitude.  Therefore,  the  effect  of  velocity  on  6e  had  the 
greatest  impact  on  the  straight  and  level  flight  results.  In  Regime  2  the  flight  test  6e  was 
consistently  0.75°  higher  than  the  predicted  curve.  The  flight  test  data  were  tightly 
grouped  around  the  average  trend-line  and  showed  the  strong  effect  that  a  change  in 
velocity  had  in  a  change  in  the  5e  value.  The  larger  difference  between  predicted  and 
flight  test  data  in  Regime  2  was  similar  to  the  results  as  discussed  above  for  5(CG). 

The  a  and  0  results  are  plotted  in  Figures  D5  and  D6.  The  flight  test  results 
matched  the  theoretical  predictions  within  ±1.0°,  except  at  the  slowest  airspeeds  in 
Regime  2.  The  data  were  in  a  nearly  continuous  curve  across  both  regimes.  There  was  a 
slight  jump  in  the  theoretical  curve  between  Regime  1  and  Regime  2  of  less  than  0.2°. 
Again  there  was  more  scatter  in  the  a  flight  test  results  than  in  the  0  values.  The  a  and  0 
flight  test  data  matched  the  model  very  closely.  The  velocity  effect  dominated  and  the 
data  were  scattered  by  the  other  factors. 

Overall,  there  was  a  strong  correlation  between  changes  in  true  airspeed  and  5e. 
The  correlation  was  also  seen  in  a  and  0.  Velocity  measurement  accuracy  was  one  of  the 
key  indicators  of  6e  prediction  accuracy.  Example  calculations  showing  the  accuracy 
using  actual  flight  conditions  are  shown  in  Appendix  E. 

Compare  6(PA)  to  model  predictions. 

A  trim  shot  was  flown  to  establish  the  conditions  outlined  for  each  point  in  the 
Test  Matrix  (see  Appendix  C).  The  test  points  were  distributed  approximately  on  the 
same  dynamic  pressure  line  (constant  airspeed),  and  have  different  Mach  numbers. 

The  altitude  effect  on  the  flight  test  parameters  was  compared  to  the  predictions. 
The  results  were  very  scattered  with  little  correlation.  Therefore  the  data  were  further 
separated  into  three  aircraft  velocity  bands.  The  flight  test  results  in  Regime  1  closely 
matched  the  predicted  values,  while  the  Regime  2  flight  test  data  were  biased  from  the 
predicted  curves. 

Figure  D7  shows  the  6e  flight  test  data  in  Regime  1  with  an  overlay  of  the 
predicted  5e  for  three  true  airspeed  values.  The  low  speed  region  (Vt  =  460  ft/s)  described 
a  much  more  shallow  curve  than  theory  predicted.  The  mid-speed  region  (Vt=  475  ft/s), 
matched  the  predicted  curve  with  a  40.2°  bias  for  all  altitudes.  The  high  speed  region  (Vt 
=  490  ft/s),  matched  the  predicted  curve,  with  about  the  same  4-0.2°  bias,  but  was  more 
scattered  than  the  mid-speed  data.  Overall  the  flight  test  data  matched  the  predicted  data 
with  a  small  amount  of  variation. 

In  Regime  2  the  flight  test  data  differed  from  the  predicted  curve  by  almost  4-1.0° 
on  average  for  all  three  velocity  bands  (see  Figure  D8).  The  flight  test  data  were  also 
more  scattered  than  in  Regime  1.  Overall  the  data  in  Regime  2  shows  the  same  general 
trend  as  the  predicted  curves,  but  the  correlation  was  not  strong.  The  altitude  effect  on  5e 
was  very  small,  and  was  masked  by  other  factors,  especially  small  changes  in  airspeed. 
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The  altitude  effect  on  a  and  0  was  also  not  very  strong.  Figures  D9  and  DIO 
show  the  plot  of  flight  test  data  in  Regime  1  for  a  and  6,  respectively,  with  the  model 
EOM  predictions  shown  for  the  same  three  true  airspeed  bands.  The  data  in  all  three 
bands  were  characterized  by  a  lot  of  scatter  in  the  flight  test  data.  All  data  were  within  2° 
of  the  predicted  curves.  As  was  discussed  in  previous  analysis,  the  6  values  were  closer 
to  the  predicted  curves  and  did  not  have  as  much  scatter  as  a.  In  Regime  2,  the  data  were 
again  scattered,  increasing  in  error  as  the  altitude  increased  (see  Figures  Dll  and  D12). 
The  flight  test  data  were  also  biased  -2°  to  -3°  from  the  predicted  value  curves.  Also  the 
a  values  were  more  scattered  than  the  0  values  at  each  altitude. 

The  altitude  effect  on  the  5‘s  was  very  small.  The  data  points  were  scattered  by 
other  effects  and  the  predicted  trends  were  difficult  to  match  to  the  model  predictions. 
The  effect  of  altitude  variations  does  not  have  a  significant  impact  on  the  result  errors. 

Compare  6(Asymmetric  T)  to  model  prediction. 

Thrust  from  one  of  the  two  engines  was  progressively  reduced  and  thrust  on  the 
other  engine  was  progressively  increased  while  keeping  the  wings  level  and  the  airspeed 
within  the  required  tolerance.  Data  were  collected  at  each  asymmetric  thrust  setting.  The 
data  were  plotted  as  the  percentage  of  total  thrust  from  the  right  engine. 

Limited  data  were  collected,  especially  in  Regime  2.  The  data  were  very 
scattered,  and  did  not  match  the  predicted  values  very  closely  for  lateral-directional  5’s 
values.  Figures  D13  and  D14  show  the  5r  and  6a  values,  respectively,  from  flight  test 
compared  to  the  predicted  theoretical  values  for  both  regimes.  The  theory  predicted  a 
larger  asymmetric  thrust  influence  on  6r  (slope  -0.5°  per  10%  change  in  asymmetric 
thrust,  verses  less  than  0.1°  per  10%  for  the  averaged  values  of  the  flight  test  data  in 
Regime  1).  The  6r  results  were  about  the  same  in  Regime  2.  The  5a  results  showed  the 
opposite  trend  result.  The  theory  predicted  about  a  -0.1°  change  per  10%  change  in 
asymmetric  T,  while  the  averaged  flight  test  data  had  a  slope  of  about  -0.5°  per  10% 
change  in  asymmetric  T.  The  data  had  a  lot  of  noise  and  stable  points  were  difficult  to 
achieve  in  asymmetrical  conditions.  The  <|)  measurements  also  had  a  lot  of  noise  and 
were  accurate  to  within  only  a  degree  or  two  of  the  averaged  value  at  most  test  points. 
These  errors  may  account  for  the  inaccuracies  between  the  flight  test  results  and  the 
theoretical  predictions. 

The  6e  values  were  tabulated  and  compared  to  the  theoretical  predictions. 
Theoretically  the  6e  was  constant  for  all  asymmetric  thrust  test  points.  Figure  D15  shows 
that  in  Regime  1 ,  all  5e  values  were  within  ±0.3°  of  the  theoretical  value.  If  the  four  data 
points  with  large  changes  in  (|)  were  discarded,  the  average  flight  test  data  6e  was  -0.03° 
from  the  theoretical  6e  and  all  flight  test  points  were  within  ±1.0°  of  the  predicted  value. 
In  Regime  2  the  flight  data  6e  value  was  biased  -1.0°  with  very  little  scatter  in  the  data. 

The  lateral-directional  5  errors  also  were  increased  because  the  test  points  were 
not  exactly  in  perfect  trim  condition  in  most  cases.  Figure  D16  and  D17  show  the  plot  of 
P  and  (j),  respectively,  for  each  test  point.  The  theoretical  values  for  the  <]>  were  expected 
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to  be  0.0°  because  the  maneuver  was  designed  for  wings  level  flight.  Flight  test  points 
were  not  perfectly  trimmed  to  wings  level,  straight  flight,  so  this  plot  is  a  measure  of  the 
error  due  to  non-ideal  flight  test  conditions.  The  errors  in  (j)  were  all  less  than  ±4.0°.  The 
offsets  from  trim  would  also  contribute  to  the  errors  in  lateral-directional  6’s.  The 
(j)  errors  were  largest  for  the  points  around  50%  asymmetric  T,  which  correspond  to  the 
largest  errors  on  the  5  plots  (Figures  D16  and  D17).  Also  the  ({)  errors  tended  to  be  more 
negative  as  asymmetric  thrust  approached  0%,  and  more  positive  as  asymmetric  thrust 
approached  100%.  This  was  also  observed  in  the  5a  and  5r  values  at  the  extreme 
asymmetrical  thrust  conditions. 

The  lateral-directional  stability  derivatives  had  a  much  larger  uncertainty  than  the 
longitudinal  stability  derivatives.  This  lateral-directional  stability  derivative  uncertainty 
contributed  to  variations  in  model  5’s  as  indicated  by  the  sensitivity  analysis.  The  most 
critical  stability  derivatives  (CnSr,  Cyp,  Cnp  and  Cip)  created  a  large  error  for  a  small 
variation.  Finally,  a  bias  was  observed  in  both  the  5a  and  5r  values.  For  Flight  10,  a 
constant  +1.6°  bias  was  observed  in  all  straight  and  level  test  cases.  The  5a  had  a  small 
bias  of  about  0.3°  for  Flight  10.  The  5a  and  6r  were  manually  trimmed,  and  the  small 
error  could  not  be  detected  from  the  cockpit.  The  trim  error  contributed  directly  to  an 
error  in  the  5’s  results. 

Overall  the  asymmetric  thrust  flight  test  data  did  not  match  the  theoretical 
predictions.  The  correct  trend  was  seen  in  the  data,  but  the  flight  test  results  had  large 
differences  from  the  expected  values.  The  difficulty  in  establishing  trimmed  flight  in  an 
asymmetrical  flight  condition  and  the  inaccuracy  in  the  bank  angle  measurement 
aggravated  the  errors.  A  sample  calculation  of  the  5’s  using  exact  flight  conditions  for  an 
asymmetrical  thrust  case  can  be  found  in  Appendix  E. 

Compare  6(7)  to  model  prediction. 

A  sawtooth  climb/descent  was  flown  for  each  pitch  angle  to  establish  the 
conditions  outlined  for  each  point  in  the  Test  Matrix  (see  Appendix  C)  to  investigate  y. 

The  Y  effect  on  5e  is  shown  in  Figure  D18.  The  y  effect  on  5e  was  predicted  to  be 
very  small,  and  the  flight  test  data  matched  the  predictions.  In  Regime  1,  all  the  flight 
test  data  were  within  ±0.25°  of  the  theoretical  curve.  In  Regime  2  the  flight  test  5e  was 
biased  +0.7°  from  the  predicted  value.  The  flight  test  data  were  also  more  scattered  about 
the  average  value  in  Regime  2  than  in  Regime  1.  Overall,  the  predicted  y  correlation  to  5e 
matched  flight  test  results,  but  the  impact  of  y  on  5e  was  very  small. 

Figures  D19  and  D20  shows  the  comparison  of  flight  test  y  effect  on  a  and 
0  respectively  to  the  theoretically  predicted  relationship.  In  Regime  1,  all  flight  test  data 
for  both  CL  and  0  were  within  ±1.0°  of  the  model  EOM  prediction  for  all  y  values.  In 
Regime  2  the  flight  test  data  were  within  ±2.0°  of  the  predicted  value,  with  a  slight  bias 
of  about  -1.0°.  As  seen  for  other  longitudinal  investigations,  the  scatter  in  the  a  data 
points  was  larger  than  the  scatter  in  the  0  values. 
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Overall,  the  flight  test  data  correlated  with  the  predicted  values  for  all  y.  The 
impact  of  y  was  not  significant,  but  was  predicted  properly  by  theory. 

Compare  6(p)  to  model  prediction. 

A  steady  heading  sideslip  was  flown  to  establish  the  conditions  outlined  for  each 
point  in  the  Test  Matrix  (see  Appendix  C).  The  maneuver  was  not  a  dynamic  maneuver. 
At  each  test  point  data  were  recorded  with  the  aircraft  in  a  stable  condition.  Data  were 
graphed,  comparing  degrees  of  6a  and  8r  to  the  steady-state  P  that  resulted  from  the 
deflections,  and  compared  P  to  (]). 

Regime  1  and  Regime  2, 6a  and  6r  were  compared  to  the  6a  and  6r  predicted  by  the 
model  EOM.  The  6r  is  shown  in  Figure  D21  and  6a  is  shown  in  Figure  1)22.  The  flight 
test  data  for  6a  and  6r  were  linear,  as  was  predicted  by  the  model.  In  addition,  the  slopes 
of  the  flight  test  data  had  similar  slopes  to  that  predicted  by  the  model.  When  the  slopes 
did  vary  from  the  model  predictions,  the  slope  of  the  flight  test  data  were  higher,  but  only 
by  approximately  0.1°  of  6  per  degree  of  p.  It  was  also  noted  that  at  P  =  0°  6a  was 
negative  and  6r  was  positive.  The  deflection  at  P  =  0°  had  much  more  effect  on  moving 
the  data  away  from  the  prediction  than  the  slopes  of  the  data.  The  rudder  showed  zero 
deflection  at  approximately  P  =  -2°,  and  the  ailerons  showed  no  deflection  at 
approximately  p  =  -0.5°.  It  is  unknown  if  this  was  an  instrumentation  issue,  a  slightly 
bent  aircraft,  or  some  other  phenomenon. 

The  (])  flight  test  data  were  also  graphed  for  the  steady-heading  sideslips,  as  seen 
in  Figure  D23.  For  Regime  1  and  Regime  2,  the  slope  from  the  flight  test  data  were 
approximately  2.5  times  greater  than  the  predicted  slope. 

Overall,  the  data  quality  from  flight  test  was  good  and  showed  a  linear 
relationship  between  6’s  and  p,  and  <])  to  p.  The  slopes  of  the  6’s  to  P  closely  matched 
model  prediction,  while  the  slope  of  (])  to  p  differed  from  the  model  by  a  factor  of  2.5. 

Compare  b(^)  to  model  prediction. 

Sustained  coordinated  level  turns  at  different  (j)  were  flown  at  the  conditions 
outlined  for  each  point  in  the  Test  Matrix  (see  Appendix  C). 

Figure  D24  shows  the  6r  and  Figure  D25  shows  6a  as  a  function  of  (])  for  Regime  1 
and  Regime  2.  Regime  1  6a  increased  approximately  0.13°  per  10°  of  bank. 
Additionally,  a  positive  6a  (stick  to  the  right)  was  required  to  maintain  a  right  bank.  This 
indicated  a  small  magnitude  of  positive  spiral  stability.  Zero  0  required  zero  6a,  which 
was  to  be  expected.  Regime  2  6r  increased  approximately  0.04°  per  10°  of  bank. 
However,  at  (|)=1°,  the  rudder  indicated  a  deflection  of  just  greater  than  1°.  This  indicated 
a  slight  bias  of  the  neutral  rudder  position.  Regime  2  6a  data  were  scattered  and  did  not 
show  a  definite  trend.  Regime  2  6r  data  were  also  scattered  with  no  trend  information. 
Part  of  the  lack  of  trend  information  available  for  Regime  2  was  that  there  were  a  very 
limited  amount  of  ^  data  in  Regime  2  and  a  very  small  spread  of  (|)’s  examined. 
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The  dynamic  flight  test  maneuvers  were  not  well  predicted  by  the  equation 
models.  The  Ig  model  assumption  was  violated.  Therefore,  the  model  was  not  able  to 
match  the  flight  test  results.  Including  the  load  factor  in  the  model  equations  may 
improve  the  prediction  accuracy.  Flight  test  data  from  less  dynamic  maneuvers  matched 
the  predicted  results  with  small  errors.  In  general  the  closer  the  trimmed  flight  test  point 
was  to  straight  and  level  flight  the  more  closely  the  flight  test  results  matched  the 
prediction. 

The  theoretical  predictions  did  not  match  the  flight  test  data  in  either  region.  This 
was  due  to  the  model  assumption  of  1-g  flight,  while  the  flight  test  was  performed  in 
level  turns  with  an  increasing  load  factor.  This  caused  the  model  to  break  down  and  a 
comparison  to  flight  test  data  useless. 

Compare  5(Q)  to  model  prediction. 

A  symmetric  push-over  was  flown  from  approximately  15°  nose-high  to  5°  nose 
low  at  different  g  (or  Q).  Data  were  collected  while  the  aircraft  passed  through 
approximately  level  with  the  horizon,  at  the  required  load  factor. 

Flight  test  data  showed  a  linear  relationship  between  5e  and  pitch  rate.  Figure 
D27  shows  that  a  trailing-edge  down  5e  resulted  in  a  pitch  over,  as  would  be  expected, 
and  a  greater  control  surface  input  resulted  in  a  higher  Q.  The  model  for  5e  showed  that 
5e  was  nearly  invariant  with  respect  to  Q.  This  was  due  to  the  fact  that  the  model  was 
written  for  1.0  g  flight.  This  assumption  makes  the  data  from  the  model  not  comparable 
to  the  flight  data  and  no  conclusions  may  be  drawn. 

The  other  model  comparison  that  the  test  team  compared  flight  data  to  and 
graphed  was  Q  and  a.  Figure  D28  shows  that  there  was  a  linear  relationship  between 
a  and  Q.  When  the  theoretical  data  were  calculated  from  the  model,  a  was  constant  for  a 
given  pitch  rate.  This,  again,  was  a  problem  with  the  model,  for  the  model  was  set  up  for 
1.0  g  flight,  where  a  given  a  would  result  in  Q=0°/sec  (level  flight).  Consequently, 
comparison  of  the  data  to  the  model  was  not  useful  without  having  the  model  incorporate 
flight  conditions  at  other  than  1.0  g.  Apply  load  factor  effects  to  the  model  EOM  and 
re-compare  to  flight  test  data  (R6). 
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Determine  the  influence  on  the  model  accuracy  of  the  non-linear  aircraft 
response  when  flying  close  to  the  stall  condition  and  the  transonic  region. 

For  each  flight  regime,  the  average  error  of  the  6  as  a  function  of  each  variable 
(W,  h,  V,  Y,  P>  Q>  asymmetric  T)  was  calculated.  Table  7  summarizes  the  results. 


Regime  1 

Regime  2 

8.  (deg) 

5e(deg) 

5a  (deg) 

5r(deg) 

h 

n/a 

n/a 

3.50 

n/a 

n/a 

0.18 

n/a 

n/a 

1.01 

n/a 

n/a 

1.24 

n/a 

n/a 

2.04 

n/a 

n/a 

AT 

0.53 

0.97 

0.99 

0.27 

phi 

1.97 

0.43 

0.21 

4.30 

1.85 

■m 

Q 

5.06 

n/a 

n/a 

2.90 

n/a 

n/a 

beta 

0.2 

0.64 

0.39 

0.23 

1.33 

gamma 

n/a 

n/a 

n/a 

0.6 

n/a 

n/a 

Table  7:  Regimes  Comparison 

The  data  showed,  for  each  regime  and  for  each  variable,  the  absolute  value  for  the 
average  difference  between  the  calculated  and  the  predicted  6’s.  Overall,  since  the  two 
regimes  were  characterized  by  a  different  range  of  a  (so  the  hypothesis  of  linear 
aerodynamic  flow  was  not  applicable  in  Regime  2),  the  model  prediction  showed  a  better 
estimation  in  Regime  1  than  in  Regime  2.  However,  only  in  a  few  cases  was  the 
difference  larger  than  1°;  in  both  cases  the  load  factor  was  different  from  1.0  g.  The 
model  equations  were  already  deemed  to  be  inapplicable  with  load  factor  different  from 
1.0  g,  as  discussed  above.  Therefore,  if  the  comparison  between  the  two  regimes  ignored 
the  two  cases  where  the  equations  were  not  applicable,  the  model  showed  a  relatively 
good  flexibility,  despite  the  different  flight  regimes. 
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CONCLUSIONS  AND  RECOMMENDATIONS 

Determination  of  Model  Sensitivity  to  Stability  Derivative  Variation 

In  the  longitudinal  axis,  Cmo,  C^a  and  Cut  were  the  most  critical  derivatives 
because  a  small  percentage  difference  in  their  estimation  yielded  a  difference  in  the 
elevator  defection  greater  than  one  degree.  This  issue  did  not  affect  the  quality  of  the 
model  predictions  because  sensitivity  analysis  performed  concluded  the  longitudinal 
derivatives  estimation  was  satisfactory. 

In  the  lateral-directional  axis,  Cjp,  Cyp,  Cisa,  CnSr,  Cnp  were  found  to  be  the  most 
critical.  More  specifically  the  rudder  control  power  was  the  most  critical  in  the  rudder 
deflection  prediction,  while  the  dihedral  effect  and  the  aileron  control  power  were  the 
most  critical  in  the  aileron  deflection  prediction.  Cyp  was  found  to  significantly  impact 
the  prediction  in  both  aileron  and  rudder  deflections.  These  results  affected  the  quality  of 
the  model  predictions  because  the  estimations  of  Cip  and  Cisa  were  marginal.  On  the 
other  hand,  Cyp  and  the  directional  derivatives  estimations  were  found  satisfactory. 

Determination  of  Test  Aircraft  Stability  Derivatives 

Overall,  the  quality  of  the  data  collected  was  satisfactory  in  the  longitudinal  and 
the  directional  axis,  and  marginal  in  the  lateral  axis.  The  variations  of  the  derivatives 
with  angle  of  attack,  sideslip,  Mach  number  and  dynamic  pressure  were  deemed 
predictable,  except  for  two  cases:  the  aileron  control  power  as  function  of  Mach  in 
Regimel,  and  the  dihedral  effect  as  function  of  dynamic  pressure  in  Regime  2. 

For  the  first  case,  it  was  possible  that  the  data  base  PEST  produced  assigned  too 
few  significant  digits  (only  six)  to  the  derivatives  and  the  lOCRP’s  which  skewed  the 
data.  For  the  second  case  no  valid  explanations  were  found. 

The  time  history  comparison  showed  that  the  angle  of  attack  and  the  angle  of 
sideslip  gauges  needed  better  compensation  for  the  time  delay  for  the  aircraft  normal, 
longitudinal  and  lateral  accelerations.  This  characteristic  did  not  significantly  affect  the 
accuracy  of  the  predictions.  However,  for  future  testing,  to  further  improve  the 
aerodynamic  model,  it  was  recommended  to  compensate  for  these  factors. 

R1  Determine  the  source  of  and  correct  for  the  discrepancy  between  PEST 
predictions  of  a  and  p  and  onboard  data  acquisition  values  of  a  and  p  prior  to 
future  stability  derivative  determination  on  the  test  aircraft.  (Page  11) 

The  zero  lift,  the  adverse  yaw  and  the  roll  due  to  the  rudder  deflection  coefficients 
were  characterized  by  significant  scatter  and  were  not  plotted.  The  lack  of  accuracy  in 
these  derivatives  was  predictable  and  further  data  analysis  was  suggested. 
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R2  Determine  the  source  of,  and  correct  for,  the  signiHcant  data  scatter  of 
Cdo,  CnSa  and  Ci5r  prior  to  future  attempts  at  determination  of  these  specific  stability 
derivatives.  (Page  12) 

The  flight  regime  seemed  to  affect  the  derivatives  mostly  because  the  angles  of 
attack  at  which  the  flight  test  techniques  were  flown  in  the  two  regimes  were  significantly 
different.  If  the  maneuvers  were  performed  at  the  same  angles  of  attack,  the  influence  of 
the  flight  regime  probably  would  have  been  negligible.  For  future  testing,  it  was 
suggested  to  normalize  the  data  collected  for  the  CG  position  and  for  the  angle  of  attack. 
Moreover,  it  was  suggested  to  fly  the  flight  test  techniques  where  Mach  number  or 
dynamic  pressure  effects  are  evaluated,  at  the  same  angle  of  attack.  For  example 
performing  the  maneuvers  in  a  wind  up  turn  could  better  isolate  Mach,  dynamic  pressure 
and  angle  of  sideslip  from  angle  of  attack. 

The  stabilator  position  was  also  changed  between  Regime  1  and  Regime  2.  The 
position  of  the  stabilator  had  an  effect  on  the  elevator  trim  values,  and  may  have  changed 
the  derivatives  between  the  two  regimes. 

R5  Investigate  the  efi'ect  of  test  aircraft  stabilator  position  on  the  stability 
derivatives.  (Page  15) 

The  Ci6a  tends  to  increase  with  M.  Theory  does  not  support  this  result;  therefore 
further  data  analysis  was  suggested.  One  possible  explanation  could  be  the  limited 
number  of  significant  digits  assigned  by  PEST  to  the  derivatives. 

Refine  Cisa  as  a  function  of  M  in  Regime  1  during  future  test  aircraft  stability 
derivative  determination  (R4). 

In  Regime  2  Cib  decreased  with  q.  The  dynamic  pressure  was  not  expected  to 
have  any  influence  on  the  stability  derivatives,  in  Regime  2. 

Refine  Cip  as  a  function  of  q  in  Regime  2  during  future  test  aircraft  stability 
derivative  determination  (R5). 

Comparison  of  Flight  Test  5’s  to  Predicted  5’s 

Overall,  the  measured  control  surface  deflections,  satisfactorily  compares  to  the 
model  equations  of  motion  predicted  surface  deflections.  The  longitudinal  cases  matched 
for  all  of  the  trim  conditions  except  for  5(Q).  The  lateral-directional  cases  did  not  show 
as  much  correlation  between  the  flight  test  data  and  the  predicted  results  as  the 
longitudinal  cases.  For  a  specific  case  with  the  exact  flight  test  conditions  input  into  the 
model  equations,  the  predicted  control  surface  deflections  matched  the  flight  test 
measurements  to  within  1.0°  for  a  majority  of  the  cases.  The  model  predictions  were  also 
more  accurate  for  straight  and  level  flight  conditions. 
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The  test  point  set-up  had  a  major  impact  on  variations  between  the  flight  test  data 
and  the  theoretical  predictions.  The  flight  test  velocity  had  the  largest  impact  on  changes 
in  the  elevator  defection  values.  Small  variations  in  test  point  velocity  produced  a  large 
change  in  elevator  deflection  values.  The  bank  angle  varied  up  to  ±8.0°  from  the 
specified  test  conditions.  Often  there  was  a  small  roll  rate  (±0.5°/sec),  that  was  not 
observed  in  flight,  but  that  impacted  the  flight  test  data.  The  bank  angle  measurement 
also  had  the  largest  oscillation  about  the  steady-state  value,  with  up  to  a  ±1.0°  oscillation 
about  the  steady-state  value.  The  combined  errors  in  bank  angle  values  had  a  significant 
impact  on  the  lateral-directional  results.  Combining  variations  resulted  in  a  combined 
change  in  the  control  surface  value,  which  accounts  for  some  of  the  scatter  seen  in  the 
flight  test  data. 

The  flight  test  data  and  equations  of  motion  predictions  matched  much  better  in 
the  longitudinal  cases  than  in  the  lateral-directional  cases.  The  variations  in  lateral- 
directional  cases  were  caused  by  the  small  variations  in  test  point  set-up,  the  higher 
uncertainty  of  the  lateral-directional  stability  derivatives,  and  the  biases  in  trim  of  the 
lateral-directional  control  surfaces.  The  lateral -directional  test  points  were  generally 
more  difficult  to  establish  in  a  steady  trim  condition  than  the  longitudinal  test  cases. 
Therefore,  the  lateral-directional  flight  test  data  tended  to  have  more  variations  from  the 
stable  trim  condition.  The  bank  angle  variations  also  caused  a  significant  impact  on 
lateral-directional  control  surface  deflection  values.  Finally,  a  bias  was  observed  in  both 
the  aileron  and  rudder  deflection  values.  For  Flight  10,  a  constant  -t-1.6°  bias  was 
observed  in  all  straight  and  level  test  cases  where  the  actual  flight  test  results  were  input. 
The  aileron  had  a  small  bias  of  about  0.3°  for  Flight  10.  The  rudder  and  aileron  were 
manually  trimmed,  and  the  small  error  could  not  be  detected  from  the  cockpit.  The  trim 
error  contributed  directly  to  an  error  in  the  control  surface  deflection  results. 

The  dynamic  flight  test  maneuvers  were  not  well  predicted  by  the  equation 
models.  The  model  assumptions  were  violated  both  bank  angle  variations  and  symmetric 
push  overs  to  generate  pitch  rate.  Therefore,  the  model  was  unable  to  match  the  flight 
test  results.  The  less  dynamic  maneuvers  flight  test  data  matched  the  predicted  results 
with  small  errors.  In  general,  the  closer  the  trimmed  flight  test  point  was  to  straight  and 
level  flight  the  more  closely  the  flight  test  results  matched  the  predictions. 

R6  Apply  load  factor  effects  to  the  model  EOM  and  re-compare  to  flight  test 
data.  (Page  24) 
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Project  recommendations  are  repeated  below  in  a  prioritized  order: 

-  Apply  load  factor  effects  to  the  model  EOM  and  re-compare  to  flight  test 
data  (R6). 

-  Determine  the  source  of  and  correct  for  the  discrepancy  between  PEST 
predictions  of  a  and  p  and  onboard  data  acquisition  values  of  a  and  p 
prior  to  future  stability  derivative  determination  on  the  test  aircraft  (Rl). 

-  Determine  the  source  of,  and  correct  for,  the  signiflcant  data  scatter  of 
Cdoj  Cn5a  and  Ci5r  prior  to  future  attempts  at  determination  of  these 
specific  stability  derivatives  (R2). 

-  Investigate  the  effect  of  test  aircraft  stabilator  position  on  the  stability 
derivatives  (R3). 

-  Reflne  Cisa  as  a  function  of  M  in  Regime  1  during  future  test  aircraft 
stability  derivative  determination  (R4). 

-  Refine  Cip  as  a  function  of  q  in  Regime  2  during  future  test  aircraft 
stability  derivative  determination  (R5). 
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APPENDIX  A  -  FLIGHT  ENVELOPE  AND  FLIGHT  LIMITS 


Figure  Al:  Flight  Envelope  and  Limits 

Figure  Al  shows  the  flight  envelope  that  was  used  during  flight  test.  The  left  line 
was  the  stall  speed  line.  The  flight  manual  indicated  that  at  13,500  lbs  gross  weight,  in 
the  clean  configuration,  108  KIAS  was  stall  speed.  For  weights  different  than  13,500  lbs, 
the  following  formula  was  used: 


Vs  =  sqrt  (0.97*nz*W), 

where  W  was  the  aircraft  gross  weight  in  lbs,  nz  was  the  load  factor  normal  to  the  flight 
path  and  the  result  was  in  knots. 

The  top  line  of  Figure  Al  indicates  the  22,000  ft  PA  line.  The  stall  warnings  for 
the  Learjet  changed  above  22,000  ft  PA,  so  in  order  to  keep  the  same  gains  and  settings 
in  the  flight  control  system,  all  FTT’s  were  flown  at  the  maximum  sortie  altitude  of 
22,000  ft.  The  minimum  altitude  was  5,000  ft  above  ground  level  (AGL)  for  any  FTT. 

The  right  line  of  Figure  Al  indicates  the  maximum  airspeed,  which  was  325 
KIAS  with  the  VSS  engaged.  The  maximum  Mach  number  with  the  VSS  on  was  0.78  M, 
which  would  not  be  reached  until  25,000  ft.  The  maximum  velocity  to  be  flown  was  310 
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KIAS,  approximately  95%  of  the  maximum  airspeed  limit.  The  maximum  airspeed  with 
the  VSS  disengaged  was  the  lesser  of  361  KIAS  or  0.78  M. 

The  load  factor  limits  was  +0.15  to  +2.8  (with  the  flight  control  system  engaged). 
The  load  factor  limits  with  the  VSS  disengaged  were  -l.Og  to  +4.0g. 

The  maximum  sideslip  angle  was  a  function  of  the  airspeed  and  the  rudder 
deflection: 


(P  limit)  =  -0.25  *  (5r  limit  for  P=0)+  (P  limit  for  6r=0) 


Table  3  defines  the  sideslip  limit  at  no  rudder  and  the  rudder  limit  with  no  sideslip 
at  different  velocities.  This  was  used  as  input  for  the  above  equation  to  calculate  the 
maximum  allowable  sideslip: 


KIAS 

P  limit  for  5r=0 

5r  limit  for  P=0 

<214 

O 

P 

30° 

225 

13.4° 

26.8° 

250 

11.0° 

22.0° 

275 

9.4° 

o 

00 

00 

300 

0 

00 

14.6° 

350 

5.6° 

10.6° 

Table  Al:  LEARJET-25  Sideslip  Limits 
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APPENDIX  B  -  FLIGHT  REGIMES 


Figure  Bl:  Flight  Regimes 


Figure  Bl  shows  the  flight  envelope  (solid  line)  and  the  flight  regimes.  The  flight 
envelope  was  divided  into  three  subsets,  each  of  them  characterized  by  a  different 
combination  of  dynamic  pressure  and  Mach  number.  The  purpose  of  dividing  the  flight 
envelope  into  three  regimes  was  to  be  able  to  determine  the  accuracy  of  the  model  in 
predicting  the  5’s  in  flight  conditions  where  the  assumption  of  linear  flow  was  no  longer 
accurate.  From  an  aerod)mamic  model  determination  point  of  view,  the  purpose  was  to 
investigate  the  combined  influence  of  q  and  M  together,  as  well  as  the  influence  of  a 
different  flow  around  the  elevator,  especially  in  the  surface  control  power  and  in  the 
damping. 

The  flight  regimes  were  defined  as: 

2)  Low  M  (0.21  -  0.38),  low  q  (0.31  psi  -  0.76  psi). 

1)  Medium  M  (0.29  -  0.57),  medium  q  (0.76  psi  -  1.53  psi). 

3)  High  M  (>0.44),  high  q  (>1.53  psi). 

Constant  q  lines  delineate  the  three  regions.  The  square  markers  indicate  the  center  flight 
condition  for  each  regime. 
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APPENDIX  C  -  TEST  CONDITION  MATRIX 


Figure  Cl;  Test  Condition  Matrix 


Figure  Cl  shows  the  entire  series  of  desired  test  points  to  be  flown.  Due  to  time 
constraints,  only  data  in  Regime  1  and  Regime  2  were  collected  during  flight  test.  The 
dotted  straight  lines  define  the  three  flight  regimes  and  replace  the  constant  dynamic 
pressure  lines  for  simplicity.  The  thick  solid  line  indicates  the  flight  envelope  of  the  test 
aircraft.  The  points  shown  cover  both  the  aerodynamic  model  determination  and  the 
equation  of  motion  validation  objective.  For  each  flight  regime  the  center  points  (square 
markers),  the  q  line  across  the  center  points  and  the  maximum  and  minimum  q  tested  are 
plotted.  The  points  along  the  center  points  q  lines  were  used  to  evaluate  the  influence  of 
Mach  number  on  the  derivatives  (constant  dynamic  pressure)  and  the  influence  of  altitude 
on  the  trimmed  5’s  (the  points  along  the  constant  center  points  indicated  airspeed  lines 
are  almost  the  same  as  the  points  along  the  constant  q  lines  because,  within  the  altitude 
range  available,  the  density  doesn’t  change  significantly;  it  was  assumed  that  they 
coincide  within  the  tolerances  fixed  for  each  sortie).  The  points  along  straight  horizontal 
lines  investigated  the  effect  of  true  airspeed  on  trimmed  deflections.  The  points  along  a 
vertical  line  investigated  the  effect  of  q  on  the  derivatives  (constant  M). 

The  square  points  indicate  the  center  test  points.  The  circles  indicate  the  points 
where  h,  V,  q  and  M  influence  were  evaluated.  The  most  important  difference  between 
these  points  is  that  the  circles  indicate  that  only  a  set  of  pitch  and  rudder/aileron  doublets 
in  SLUF  were  required.  The  5’s  were  estimated  as  a  function  of  h  and  V,  while  the 
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aerodynamic  stability  derivatives  were  estimated  as  a  function  of  q  and  M.  At  the  square 
points  the  block  of  maneuvers  to  be  flown  was  more  complex.  The  pitch  and 
rudder/aileron  doublets  were  applied  in  a  SLUF  trim  condition,  in  a  climb/descent,  in  a 
stable  turn  at  different  load  factors,  in  a  push-over,  in  an  asymmetric  thrust  investigation, 
and  in  a  stable  condition  with  steady  sideslip. 

See  the  Test  Matrix  in  Table  Cl  for  more  details. 
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25  Elevator  Deflection  Velocity  Effect,  Project  HAVE  TRIM 
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Elevator  Deflection  Altitude  Effect,  Project  HAVE  TRIM 
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ieure  D8:  Learjet  25  Regime  2  Elevator  Deflection  Altitude  Effect,  Project  HAVE  TRIM 
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Figure  D9:  Learjet  25  Regime  1  Angle  of  Attack  Altitude  Effect,  Project  HAVE  TRIM 
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Figure  D12:  Learjet  25  Regime  2  Pitch  Angle  Altitude  Effect,  Project  HAVE  TRIM 
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Figure  D15:  Learjet  25  Elevator  Deflection  Asymmetric  Thrust  Effect,  Project  HAVE  TRIM 
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7:  Learjet  25  Bank  Angle  Asymmetric  Thrust  Effect,  Project  HAVE  TRIM 
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Elevator  Deflection  Flight  Path  Angle  Effect,  Project  HAVE  TRIM 
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Attack  Flight  Path  Angle  Effect,  Project  HAVE  TRIM 
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;ure  D20:  Leaijet  25  Pitch  Angle  Flight  Path  Angle  Effect,  Project  HAVE  TRIM 
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25  Bank  Angle  Sideslip  Effect,  Project  HAVE  TRIM 
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•e  D24:  Learjet  25  Rudder  Deflection  Bank  Angle  Effect,  Project  HAVE  TRIM 
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Figure  D27:  Learjet  25  Elevator  Deflection  Pitch  Rate  Effect,  Project  HAVE  TRIM 
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APPENDIX  E  -  SAMPLE  MODEL  PREDICTIONS 


General 

The  model  predictions  were  calculated  for  general  test  conditions  to  allow  a 
comparison  between  flight  test  data  and  the  model  predictions.  Small  variations  in  the 
test  point  set-up  caused  significant  impacts  in  the  variation  between  flight  test  results  and 
the  model  predicted  values.  All  test  cases  used  data  that  had  as  little  noise  as  possible 
and  flight  test  points  with  the  smallest  variations  from  the  stable  value.  In  Capt  Miller’s 
Master’s  thesis,  individual  test  point  conditions  were  investigated  to  determine  the 
accuracy  between  the  flight  test  data  and  the  model  using  the  actual  flight  test  data.  The 
point-by-point  investigation  determined  the  accuracy  that  the  model  can  predict  for  the 
required  control  surface  deflection.  Example  calculations  are  included  below  for 
comparison.  At  each  test  point,  the  calculated  stability  derivatives  and  the  actual  aircraft 
weight,  altitude,  airspeed,  flight  path  angle,  sideslip  angle  asymmetric  thrust  level  and 
moments  of  inertia  were  input  into  the  model.  The  model  predictions  for  aircraft  state 
(angle  of  attack,  pitch  angle,  thrust  settings,  bank  angle  and  angular  rates)  were  calculated 
and  then  the  control  surface  deflections  were  determined.  The  predictions  were  then 
compared  to  the  flight  test  results. 

Test  Point  87  (Flight  10):  The  test  point  was  at  straight  and  level  flight,  flown  at  the 
center  test  position  (15,000  feet,  velocity  475  feet/second)  with  a  reduced  aircraft  weight 
of  11,775  pounds.  The  predicted  aircraft  states  were  close  to  the  expected  measured 
flight  test  values.  Angle  of  attack  was  predicted  to  be  3.900  degrees,  which  was  1.03 
degrees  greater  than  the  flight  test  data.  The  predicted  pitch  angle  was  4.002  degrees, 
which  was  0.66  degrees  higher  than  the  measured  value.  The  predicted  bank  angle  was  - 
0.016  degrees,  0.13  degrees  higher  than  the  measured  flight  test  result.  The  thrust 
settings  were  both  within  a  pound  for  both  engines.  The  predicted  angular  rates  were 
within  0.09  degrees  per  second  for  all  three  axes,  with  the  yaw  rate  prediction  exactly 
matching  the  flight  test  data.  The  elevator  deflection  prediction  was  -0.507  degrees, 
which  was  only  0.11  degrees  greater  than  the  flight  test  data  and  about  equal  to  the 
measurement  precision  of  0.1  degree.  The  aileron  deflection  prediction  (0.362  degrees) 
was  within  0.03  degrees  of  the  flight  test  value.  The  rudder  deflection  prediction  was 
-0.383  degrees,  which  was  1.62  degrees  higher  than  the  flight  test  results.  A  1.6  degree 
bias  was  observed  in  all  straight  and  level  test  points  for  Flight  10,  therefore  the  bias  was 
probably  caused  by  an  asymmetrically  trimmed  rudder. 

Test  Point  60  (Flight  10):  The  test  point  was  at  a  trim  condition  for  the  bank  angle 
investigation,  and  was  flown  at  the  center  test  position  (15,000  feet,  475  feet  per  second 
at  12,377  pounds)  and  with  zero  bank  angle.  Again  the  predicted  aircraft  states  closely 
matched  the  flight  test  data.  The  predicted  angle  of  attack  was  4.04  degrees,  which  was 
0.51  degrees  below  the  flight  test  value.  The  predicted  pitch  angle  was  3.959  degrees, 
which  was  0.38  degrees  higher  than  the  measured  results.  The  model  bank  angle  was 
-0.248  degrees,  which  was  within  +0.14  degrees  of  the  flight  test  data.  The  predicted 
flight  test  thrust  was  slightly  less  than  the  measured  values,  with  just  over  10  pounds 
more  thrust  on  each  engine.  As  in  the  previous  case,  the  angular  rates  were  within  0.07 
degrees  per  second  for  all  three  axes,  with  the  raw  rate  prediction  exactly  matching  the 
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flight  test  data.  The  elevator  deflection  prediction  was  -0.658  degrees,  only  0.17  degrees 
lower  than  the  measured  flight  test  results.  The  aileron  deflection  prediction  of  0.618 
degrees  was  0.34  degrees  higher  than  the  flight  test  value.  Finally  the  rudder  deflection 
prediction  was  -0.887  degrees,  which  was  1.69  degrees  lower  than  the  flight  test  results. 
The  bias  was  nearly  equal  to  the  bias  in  Test  Point  87  of  Flight  10. 

Test  Point  27  (Flight  10):  The  test  point  corresponds  to  a  10%  asymmetric  thrust  point  at 
the  center  test  position.  The  point  was  a  lateral-directional  test  point  that  incorporated 
more  maneuvering  flight  than  the  straight  and  level  stable  condition.  The  aircraft  state 
results  matched  closely,  while  the  control  surface  deflections  were  invalid.  The  predicted 
angle  of  attack  was  4.397  degrees,  0.46  degrees  lower  than  the  flight  test  data.  The 
predicted  pitch  angle  was  4.585  degrees,  which  was  0.81  degrees  above  the  measured 
value.  The  model  bank  angle  (3.552  degrees),  differed  significantly  (+2.21  degrees)  from 
the  flight  test  results.  The  predicted  thrust  was  slightly  higher  than  the  measured  values 
(about  40  pounds  on  the  left  engine  and  15  pounds  on  the  right  engine).  Also  in  this  case, 
the  aircraft  had  larger  angular  accelerations  (the  largest  was  0.275  degrees  per  second  for 
yaw  rate),  but  the  errors  were  still  under  0.09  degrees  per  second  with  a  nearly  perfect 
match  in  yaw  rate.  The  control  surface  deflection  predictions  all  exceeded  the  maximum 
deflection  limits  and  were  therefore  invalid.  Test  Point  27  was  one  of  the  worst  case 
scenarios,  and  will  be  investigated  to  attempt  to  determine  the  source  of  the  error. 
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APPENDIX  F  -  PLOTS  SUPPORTING  DETERMINATION  OF  TEST 
AIRCRAFT  STABILITY  DERIVATIVES 
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Longitudinal  Time  History  Comparison,  Project  HAVE  TRIM 
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ire  F3:  Learjet  25  Clu  vs  a,  Project  HAVE  TRIM 
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Figure  F7:  Leaijet  25  Cma  vs  Mach  Regime  1,  Project  HAVE  TRIM 
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Figure  FIO:  Learjet  25  Cta  vs  Mach  Regime  2,  Project  HAVE  TRIM 
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Figure  F15:  Learjet  25  C^se  vs  Mach  Regime  1,  Project  HAVE  TRIM 
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Figure  F19:  Leaijet  25  Cta  vs  q  Regime  1,  Project  HAVE  TRIM 
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e  F23:  Leaijet  25  Cmq  vs  q  Regime  1,  Project  HAVE  TRIM 
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Figure  F25:  Learjet  25  Cm5e  vs  q  Regime  1,  Project  HAVE  TRIM 
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Figure  F26;  Learjet  25  Cm5e  vs  q  Regime  2,  Project  HAVE  TRIM 
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Figure  F30:  Learjet  25  Cnr  vs.  |3  Regime  1,  Project  HAVE  TRIM 
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Figure  F31;  Leaijet  25  Cnr  vs.  |3  Regime  2,  Project  HAVE  TRIM 
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Figure  F32:  Learjet  25  CnSr  vs.  P  Regime  1,  Project  HAVE  TRIM 
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VS.  3  Regime  1,  Project  HAVE  TRIM 
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VS.  6  Regime  2,  Project  HAVE  TRIM 
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Figure  F40:  Learjet  25  C„p  vs.  Mach  Regime  1,  Project  HAVE  TRIM 
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Figure  F44:  Learjet  25  CnSr  vs.  Mach  Regime  1,  Project  HAVE  TRIM 
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Figure  F48;  Learjet  25  Cip  vs.  Mach  Regime  1,  Project  HAVE  TRIM 
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Table  Gl:  Longitudinal  Derivatives  -  Regime  1 
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APPENDIX  H  -  PLOTS  SUPPORTING  DETERMINATION  OF 
MODEL  SENSITIVITY  TO  STABILITY  DERIVATIVE  VARIATION 
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Figure  H2:  Leaijet  25  Cyp  vs  a,  Project  HAVE  TRIM 
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APPENDIX  I  -  THESIS  FLIGHT  CONTROL  EQUATIONS 

The  following  equations  are  from  the  unpublished  Master’s  thesis  by  Capt  G. 
Miller  and  are  used  to  predict  5’s.  They  are  collectively  referred  to  as  the  “model  EOM” 
or  simply,  the  “model”  in  the  HAVE  TRIM  Final  Technical  Information  Memorandum. 


Longitudinal  control  surface  equation: 


6e= 


C7(ZjCosaj  +  x.pSinaj)(T|  +  T2)  C„,^PR  (R  P  )  +C^^a  +  C^^d  +  C,„^Q 


n>5. 


where: 


1 


xt  =  horizontal  (x-axis)  distance  from  thrust  center  to  center  of  gravity 
zt  =  vertical  distance  from  thrust  center  to  center  of  gravity 
Ti  =  Thrust  from  right  engine 
T2  =  Thrust  from  left  engine 
q  =  Dynamic  pressure 
c  =  Wing  chord 

ttr  =  Angle  of  thrust  from  Xbody  axis 
lyy  =  Moment  of  inertia  about  ybody  axis 


Lateral-Directional  control  surface  equations: 

+  C4C,^^  ) 


qb[(c3C^  +C4C„^  +C4C,^  )  -i-C4C„^  )(C9C„^  +C4C,^  )] 


where: 

Cp  =  C,^  PQ  +  Cj^QR  +  CjyTsinajCrj  -  T, )  +  C4yTCOsaT(T2  -  Tj)  + 

+  ib[c3(C,  P  +  C,  R  +  C,^p)  +  C4(C„^P  +  C„  R  +  C„^P)] 

c,  =  C„^  PQ  -  C„  QR  +  C5  y^  cos  (T2  -  T, )  +  C4yT  sin  Ut  (Tj  -  T,  )  + 

+  ib[cXC„^P  +  C„  R  +  C„^P)  +  C4(C.^P  +  C,  R  +  C,^p)] 


XX  zz  xz 


C4  = 


I  I  -I 

XX  zz  xz 
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■^XX^ZZ  *xz 

yr  =  horizontal  (y-axis)  distance  from  thrust  center  to  center  of  gravity 
b  =  Wing  span 

Ixx  =  Moment  of  inertia  about  Xbody  axis 

Izx  =  Moment  of  inertia  about  Zbody  axis 

Ixz  =  Product  of  inertia  about  Xbody  axis  and  Zbody  axis 

Ci  pq  =  Product  derivative 

Ci_qr=  Product  derivative 

Cn_pq  =  Product  derivative 

Cn  qr=  Product  derivative 

The  following  assumptions  were  made  to  derive  the  equations: 

-  The  sum  of  the  forces  on  the  vehicle  are  zero,  so  there  is  no  translational 
acceleration. 

-  The  sum  of  the  moments  on  the  vehicle  are  zero,  so  there  is  no  angular 
acceleration. 

-  All  aerodynamic  forces  are  linear.  The  equations  do  not  account  for  aerodynamic 
effects  near  stall  and  at  high  speed,  when  aerodynamic  compressibility  effects 
occur. 

-  Longitudinal  and  lateral-directional  equations  are  decoupled. 

-  The  mass  properties  of  the  vehicle  do  not  change  during  the  maneuver.  The 
center  of  gravity  and  moments  of  inertia  remains  constant. 

-  Ideal  atmosphere:  no  gusts  or  disturbances,  and  the  atmosphere  can  be  modeled 
by  the  standard  atmosphere. 

-  The  stability  derivatives  are  constant  throughout  the  trim  point. 
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APPENDIX  J  -  PARAMETERS  LIST 

Aerodynamic  Parameters  (symbol)  [units] 
Aircraft  Weight  (W)  [lbs.] 


Resolution  Uncertainty 

*From  fuel  flow  and  zero  fuel  wt. 


Aircraft  Airspeed  (V)  [kts] 

0.0305 

±0.5% 

Pressure  Altitude  (h)  [ft] 

19.5 

±0.4%+25  ft 

Asymmetric  Thrust  (AT)  [%  right  engine] 

*From  Engine  Pressure  Ratio 

Glide  Path  Angle  (y)  [deg] 

*From  a  and  0 

Sideslip  Angle  (P)  [deg] 

0.00098 

±1.0% 

Aircraft  Turn  Rate  (co^)  [deg/s] 

*From  P,  Q  and  R 

Aircraft  Pitch  Rate  (coe)  [deg/s] 

*From  P,  Q  and  R 

Angle  of  Attack  (a)  [deg] 

0.00098 

±1.0% 

Pitch  Angle  (0)  [deg] 

0.00439 

±1.0  deg 

Bank  Angle  ((j))  [deg] 

0.00439 

±1.0  deg 

Acceleration  in  X-body  Axis  (nx)  [g’s] 

0.00005 

±0.07% 

Acceleration  in  Y-body  Axis  (ny)  [g’s] 

0.00005 

±0.07% 

Acceleration  in  Z-body  Axis  (nz)  [g’s] 

0.00024 

±0.07% 

Rotation  about  X-body  Axis  (P)  [deg/s] 

0.00488 

±2.3% 

Rotation  about  Y-body  Axis  (Q)  [deg/s] 

0.00195 

±2.3% 

Rotation  about  Z-body  Axis  (R)  [deg/s] 

0.00195 

±2.3% 

Control  Surface  Deflections  (symbol)  [units! 

Resolution 

Uncertainty 

Aileron  Surface  Deflection  (6a)  [deg] 

0.00195 

±1.0% 

Elevator  Surface  Deflection  (6e)  [deg] 

0.00098 

±1.0% 

Rudder  Surface  Deflection  (6r)  [deg] 

0.00195 

±1.0% 

Longitudinal  Stability  Derivatives  (symbol)  [units! 

Initial  Value 

Lift  coefficient  slope  (Cta)  [1/rad] 

5.386 

Zero  lift  coefficient  (Clo)  [— ] 

0.097 

Zero  lift  drag  coefficient  (Cdo)  [— ] 

0.0235 

Zero  lift  pitching  moment  coefficient  (Cmo)  [— ] 

0.031 

Pitch  stability  (Cma)  [1/rad] 

-0.7792 

Pitch  damping  (Cmq)  [sec/rad] 

-11.918 

Pitch  control  power  (CmSe)  [1/rad] 

-0.8995 
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Lateral-Directional  Stability  Derivatives  (symbol)  [units] 

Dihedral  effect  (Cip)  [17rad] 

Lateral  damping  (Cip)  [sec/rad] 

Lateral  control  power  (Ci5a)  [1/rad] 

Weather-cock  stability  (Cnp)  [1/rad] 

Directional  damping  (Cnr)  [sec/rad] 

Directional  control  power  (CnSr)  [1/rad] 

Adverse  yaw  derivative  (CnSa)  [l/rad] 

Adverse  roll  due  rudder  deflection  (Cisr)  [l/rad] 

Roll  coupling  due  to  yaw  rate  (Cir)  [l/rad] 

Yaw  coupling  due  to  roll  rate  (Cnp)  [l/rad] 

Yaw  derivative  due  to  sideslip  (Cnp)  [l/rad] 


Initial  Value 

-0.0802 

-0.4160 

-0.0596 

0.0900 

-0.1914 

-0.0487 

0.0011 

0.0195 

0.311 

-0.0395 

0.090 
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APPENDIX  K  -  WORKING  WITH  PEST 

Overall 


PEST  was  just  one  of  the  four  applications  of  the  Parameter  Identification 
software  package.  In  order  to  run  PEST  and  calculate  the  derivatives  four  steps  were 
followed:  project  manager  shield,  time  history  editor  shield,  parameter  estimation  shield 
and  output  generator  shield.  In  the  next  sections  these  steps  will  be  briefly  described  and 
the  lessons  learned  highlighted^.  Each  data  point  run  required  about  20  minutes. 

Project  Manaeer 

The  Project  Manager  allowed  for  setting  the  constants  (the  predicted  value  for  the 
derivatives)  the  instrumentation  correction  and  the  variables  from  which  the  aerodynamic 
model  was  calculated.  PEST  did  not  actually  need  a  good  prediction  for  the  derivatives, 
but  these  predictions  were  used  as  initial  values  to  start  the  iterations  to  calculate  the 
stability  derivatives.  This  iteration  process  stopped  when  a  minimum  in  the  specific  cost 
function  was  reached.  A  bad  initial  value  could  cause  PEST  to  converge  to  a  local 
minimum  of  the  cost  function  instead  of  the  actual  absolute  minimum.  Therefore,  since  a 
good  set  of  stability  derivatives  was  not  available  except  for  one  flight  condition,  it  was 
decided  to  use  the  derivatives  related  to  that  condition  as  a  prediction  for  all  the  PEST 
runs. 


The  instrument  corrections  page  allowed  entering  the  compensations  for  CG 
position  and  aircraft  accelerations.  Throughout  the  runs  it  was  observed  that  PEST  gave 
a  relatively  high  importance  to  the  instrument  corrections,  and  if  the  time  histories  did  not 
match,  it  was  sometimes  related  to  an  incorrect  compensation  for  CG  position  and 
accelerations.  The  constants  page  regarded  mass  properties  (moments  of  inertia,  total 
aircraft  mass),  airplane  constants  (chord,  wing  span)  and  CG  position.  This  last 
parameter  was  related  to  the  instrument  corrections  page.  Whatever  CG  position  the 
instrument  corrections  were  related  to,  that  position  was  constant  for  the  entire  test 
program.  Therefore  the  shift  between  the  current  position  and  the  reference  position  was 
entered  in  the  constant  page.  The  response  page  allowed  for  selecting  the  variables  the 
derivatives  were  calculated  from  and  the  weight  of  each  of  those  variables.  The  initial 
plan  was  to  run  PEST  a  few  times,  find  the  weights  for  each  run,  calculate  the  arithmetic 
averages  of  the  weights  and  take  them  as  constant  for  the  rest  of  the  runs.  However 
PEST  was  more  sensitive  with  regard  to  the  weights.  Thus,  the  weights  were  calculated 
and  fine  tuned  for  each  run  instead  of  using  averaged  values.  This  was  a  very  time 
consuming  process. 

The  instrument  corrections  were  assumed  to  be  zero:  the  LEARJET-25  instrument 
readings  were  corrected  for  the  CG  position;  the  compensation  for  the  aircraft 
accelerations  were  unknown,  so  they  were  set  to  zero. 


’  For  further  details  see  reference  2,  PCPID  Users  manual 
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•  The  response  biases  were  all  set  to  zero  and  “false”  except  the  variables  “qbias” 

and  “anbias”  for  the  longitudinal  case,  and  “pbias”  and  “aybias”  for  the  lateral- 
directional  case. 

Figures  K1  and  K2  are  taken  from  the  “response”  page;  they  show  an  initial  value 
for  the  weightings  (column  3),  used  in  the  longitudinal  and  lateral-directional  runs 
respectively.  Whenever  the  variable  was  not  active  (“false”),  the  weight  was  zero. 
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Figure  Kl:  Response  Page  (longitudinal  case) 
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Figure  K2;  Response  Page  (lateral-directional  case) 

Time  history  editor 

The  time  history  editor  allowed  for  the  selection  of  the  time  slice  from  which  the 
derivatives  were  calculated.  A  lesson  learned  for  the  test  team  was  that  PEST  truncated 
the  time  history  if  too  many  points  were  selected,  resulting  in  a  bad  estimation  of  the 
derivatives.  Therefore,  if  the  time  spacing  between  the  PTI  rudder  doublets  and  the  PTI 
aileron  doublets  was  too  long,  the  resulting  time  history  was  not  used.  The  time  spacing 
that  was  found  to  permit  a  better  estimation  of  the  derivatives  was  1  second. 


PEST 


The  most  important  checks  accomplished  before  actually  running  PEST  was  to 
make  sure  the  variables  that  PEST  was  asked  to  calculate  the  derivatives  from,  were 
active.  The  way  to  perform  the  check  was  to  read  the  initial  dialogue  lines  that  PEST 
displayed  and  insure  that  all  of  the  active  variables  would  not  show  up  in  the  list  of 
unknown  channels.  If  one  of  the  variables  selected  in  the  response  page  of  the  PSF 
manager  shield  was  on  that  list,  it  was  necessary  to  open  an  ASCII  file  called  “plot  setup” 
to  determine  the  active  variables.  Therefore,  before  running  PEST,  this  file  was  opened 
and  a  cross-check  between  the  PSF  manager  page  and  the  lines  of  this  file  was 
performed. 

Output  £enerator 

In  the  output  generator,  the  derivatives  that  were  just  calculated  were  displayed  in 
the  “text  output”  page,  together  with  the  lOCRP.  Also,  the  computed  and  flight  data  time 
histories  were  plotted.  The  characteristics  used  to  evaluate  the  accuracy  in  the  estimation 
were  the  peak-to-peak  time  spacing  and  amplitude.  The  derivatives  plot  displayed  the 
derivatives  that  were  just  calculated  with  the  95%  confidence  interval.  If  the  estimation 
was  deemed  satisfactory  then  the  entire  set  of  data  were  saved  in  a  data  base;  for  each  run 
the  flight  conditions,  the  derivatives  and  the  CR  bounds  were  saved  automatically  in  the 
data  base.  In  order  to  get  the  95%  confidence  for  each  derivative,  the  bounds  had  to  be 
multiplied  by  10.  Another  lesson  learned  was  that  the  data  base  assigned  only  a  6-digit 
space  for  each  number.  If  the  derivatives  or  the  CR  bounds  were  less  than  0.00001  the 
value  displayed  on  the  data  base  was  0.  Thus,  if  further  data  analysis  was  required,  a 
cross-check  with  the  test  output  page  for  each  run  was  necessary  as  was  a  manual  data 
entry  of  the  actual  value. 
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APPENDIX  L  -  PROGRAMMED  TEST  INPUT  DESCRIPTION 


Single  PTI 


Figure  LI:  Programmed  Test  Input 

Figure  LI  shows  an  example  of  doublet  PTI.  The  rise,  transition,  amplitude  and 
decay  slope  were  set  on  board  and  could  be  changed.  The  symmetry  of  the  PTI  in  time 
was  changed,  as  shown  above,  in  order  to  make  the  positive  part  different  from  the 
negative.  For  a  double  PTI,  the  one  used  for  lateral/directional  derivatives  investigation, 
the  time  between  the  rudder  doublet  and  the  aileron  doublet  could  be  varied  too. 
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APPENDIX  M  -  GROUND  SIM  AND  PTI  SHAPING 

The  test  team  could  have  started  flight  test  sooner  after  the  aircraft  arrived  on 
station  if  more  ground  work  had  been  done  in  advance.  The  test  team  had  coordinated 
what  variables  were  needed  to  perform  the  data  reduction  and  obtain  the  trim  conditions 
and  stability  derivatives.  After  the  aircraft  arrived,  the  test  team  discovered  that  the  exact 
settings  and  names  of  the  variables  had  to  match  precisely  for  the  data  to  transfer  between 
the  program  the  aircraft  used  to  save  data  and  PEST.  The  test  team  spent  considerable 
time  changing  variable  names  and  order,  such  that  PEST  could  recognize  the  data  files. 
A  sample  file  from  the  contractor  previous  to  the  aircraft  arriviving  would  have  been  very 
useful.  The  test  team  could  have  reviewed  the  file  and  check  that  the  proper  variables 
were  available  to  PEST  and  that  they  were  in  the  right  order.  The  test  team  took  a  full 
additional  day  to  get  the  data  issues  worked  out  before  flight  testing. 

PTI  shaping  was  another  issue  that  the  test  team  faced.  The  test  team  found  no 
way  to  predict  what  the  PTI  should  look  like  prior  to  flight,  in  terms  of  either  magnitude 
or  symmetry.  Symmetry  was  adjusted  such  that  the  aircraft  completed  the  maneuver  at  as 
close  to  the  initial  trim  condition  as  possible.  Symmetry  was  measured  in  percent,  from 
0.00  to  1.00.  Symmetry  ranged  from  .45  to  .48  for  test  maneuvers.  This  corresponded  to 
a  shorter  initial  input  followed  by  a  longer  second  input.  The  longer  second  input 
presumably  helped  overcome  the  rates  generated  by  the  first  input. 

The  following  discussion  on  PTI  amplitude  applies  most  specifically  to  the 
LEARIET-25,  but  may  have  some  practical  application  to  other  flight  testing.  The  test 
team  had  to  choose  amplitudes  that  would  get  measurable  responses,  but  not  trip  off  the 
variable  stability  system.  Problems  were  encountered  in  Regime  1  at  high  beta  when 
inputting  the  yaw  doublet  and  under  increased  G  loading  in  the  constant  altitude  turn. 
The  test  team  had  planned  to  test  up  to  9  degrees  of  sideslip,  approaching  the  VSS 
sideslip  limit.  Actual  testing  showed  that  it  was  only  practical  to  go  to  4°  sideslip  and 
amplitude  of  the  PTI  was  reduced  to  half  of  the  amplitude  of  the  other  conditions.  The 
rolling  moment  safety  trip  kicked  off  the  VSS  system  prior  to  the  maximum  sideslip. 
Test  points  at  60°  bank  were  accomplished  with  a  20%  reduction  in  the  PTI  pitch 
amplitude,  and  66°  bank  points  were  not  attempted.  Table  N1  is  a  summary  of  the 
amplitudes  and  symmetries  used.  Regime  2  did  not  need  to  have  the  doublets  adjusted. 
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Regime  1 

Pitch 

Roil 

Yaw 

SLUF 

Amp 

7 

15 

10 

Period 

0.6  sec 

0.8  sec 

1 .0  sec 

Sym 

0.45 

0.47 

0.49 

Banked  turns 

Amp 

4 

10 

7 

Period 

0.6  sec 

0.8  sec 

1 .0  sec 

Sym 

0.45 

0.48 

0.46 

Push  Over  0.4g 

Amp 

4 

10 

10 

Period 

0.6  sec 

0.8  sec 

1 .0  sec 

Sym 

0.45 

0.48 

0.46 

Regime  2 

Pitch 

Roll 

Yaw 

SLUF 

Amp 

7 

15 

10 

Period 

0.6  sec 

0.8  sec 

1 .0  sec 

Sym 

0.45 

0.45 

0.45 

Table  Ml:  PTI  Summary 


The  final  issue  addressed  by  the  test  team  was  the  lat-dir  doublets  were  initially  too  long 
to  be  analyzed  by  PEST.  The  end  of  the  maneuver  was  cut  off,  losing  important  data. 
The  solution  was  to  sequence  the  doublets  such  that  by  one  test  conductor  input,  the  roll 
doublet  was  triggered  followed  shortly  afterwards  (1  second)  by  the  yaw  doublet.  This 
was  determined  to  be  the  best  technique  for  getting  the  lat-dir  stability  derivatives. 
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APPENDIX  N  -  IN  FLIGHT  LESSONS  LEARNED 

General 

Record  the  cockpit  intercom  for  qualitative  comments  and  data. 

Any  flight  135  KIAS  or  slower  should  be  performed  at  a  maximum  of  3000  lbs  of  fuel 
remaining  to  allow  the  VSS  to  remain  engaged. 

Configuration 

Ensure  that  the  correct  stabilator  position  is  verified  after  each  engagement  of  the 

VSS. 

Ensure  that  the  stabilator  is  set  to  Manual  during  VSS  engagement  so  that  it  does  not 
automatically  trim. 

Return  stabilator  trim  to  Auto  during  administrative  phases  of  flight.  Otherwise,  the 
stabilator  will  cause  a  fairly  violent  pitch  up  or  down  after  a  VSS  trip  if  displaced  from 
the  trim  airspeed. 

Check  the  flaps  are  in  the  full  up  position  periodically  because  they  tend  to  droop. 

Only  the  two  heavier  CG  subsets  were  attainable. 

Asymmetric  Thrust  Points 

The  maximum  achievable  thrust  asymmetry  in  level  flight  was  around  30%.  Idle 
occurs  around  61%. 

Sideslip 

Use  the  safety  pilot  to  trim  the  rudder  to  establish  the  desired  sideslip  angle. 
Maximum  sideslip  was  around  5  degrees  before  the  VSS  would  disengage  due  to  side 
force. 

Begin  the  rudder  doublet  away  from  the  sideslip  to  help  ensure  the  VSS  remains  engaged. 
Decrease  the  PTI  rudder  doublet  amplitude  when  approaching  large  sideslip  values  (we 
decreased  the  magnitude  from  9  to  4.5). 

An  engine  flamed  out  during  the  recovery  from  a  prolonged  sideslip  with  500  lbs  of  fuel 
remaining  in  each  wing  tank  (1500  lbs  total  fuel).  The  Veridian  experts  decided  that  the 
flameout  was  due  to  fuel  starvation  with  less  than  full  fuel  in  the  wing  fuel  tanks.  The 
crew  limited  future  sideslips  to  greater  than  3000  lbs  of  fuel  remaining  (full  wing  tanks). 
No  subsequent  problems  were  discovered. 

Climbs  and  Descents 

Maximum  dive  angle  was  approximately  3  degrees  with  idle  thrust. 

Maximum  climb  angle  was  approximately  10  degrees  with  maximum  thrust. 
Approximately  3000  feet  of  altitude  was  required  to  establish  a  trimmed  climb. 
Approximately  1500  feet  of  altitude  was  required  to  establish  a  trimmed  idle  descent. 

Turns 

In  Region  1,  the  maximum  bank  angle  possible  prior  to  VSS  automatic 
disengagement  was  60  degrees. 

In  Region  2,  the  maximum  bank  angle  possible  prior  to  VSS  automatic  disengagement 
was  30  degrees. 
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Decrease  the  PTI  elevator  doublet  amplitude  when  approaching  large  bank  angle  values 
(we  decreased  the  magnitude  from  5  to  4). 

Pitch  Rate 

Hand  fly  pitch  rate  points. 

In  order  to  enter  the  pitch  rate  points  begin  at  level  flight  about  50  KIAS  above  the  data 
band.  Establish  a  12  degree  climb  using  the  safety  pilot  to  set  and  adjust  the  throttles  as 
necessary.  At  about  15  KIAS  above  the  data  band,  begin  the  bunt.  We  gained  about  500 
feet  from  the  initiation  of  the  bunt  to  the  apex  altitude.  Perform  the  PTI  approaching 
wings  level. 

PTI  and  PEST 

Allow  a  maximum  of  2  seconds  between  aileron  and  rudder  doublets.  The 
smaller  the  pause  between  doublets,  the  easier  the  data  is  to  analyze  in  PEST. 

Perform  3  of  each  type  of  doublet  at  each  condition  consecutively.  This  allows  for  easier 
PEST  analysis  because  weightings  do  not  have  to  be  changed  between  each  data  point. 
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APPENDIX  O  -  POINTS  OF  CONTACT 


Requesting  Agency 

Department  of  Aeronautics  &  Astronautics, 

Air  Force  Institute  of  Technology  (AFIT/ENY) 
2590  P  Street,  Suite  201 
Wright-Patterson  AFB,  OH  45433-7765 


Agency  Representative:  Lt  Col  David  Jacques,  USAF  AFIT/ENY 

DSN  785-3636  ext  4636 
E-mail:  david.iacQues@afit.af.edu 


Responsible  Test  Organization 

The  RTO  for  this  project  is  the  412**’  Test  Wing.  The  HAVE  TRIM  team  will 
accomplish  all  testing  and  evaluation.  The  test  team  consists  of  two  student  test  pilots, 
one  student  flight  test  navigator  (FTN),  and  two  student  flight  test  engineers  (FTEs). 

Project  Manager 

Captain  Gary  Miller 
USAF  TPS/EDA 

DSN  527-9968/9970/3000,  COMM  (661)  277-9968/9970/3000 
E-mail:  Gary.Miller4@edwards.af.niil 

Advanced  Information  Engineering  Services  (formerly  Veridian) 

Mr  Russ  Easter 

Advanced  Information  Engineering  Services  Instructor  Pilot  (IP) 

COMM  (716)  310-0278 

Email:  jeaster@buffalo.veridian.com 

USAF  TPS 
Mr  Russell  Erb 
USAFTPS/EDP 

DSN  527-8829,  COMM  (661)  277-8829/3000 
E-mail:  russell.erb@edwards.af.mil 


NASA  Dryden 
Mr  Chris  Nagy 

Aerospace  Engineer  /  NASA  Dryden 
COMM  (661)  276-2626 
E-mail:  @dfrc.nasa.gov 
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LIST  OF  ABBREVIATIONS,  ACRONYMS,  AND  SYMBOLS 

AFFTC  -  Air  Force  Flight  Test  Center 
AFIT  -  Air  Force  Institute  of  Technology 
AFRL  -  Air  Force  Research  Laboratory 
Alpha  (a)  -  Angle  of  Attack 
Beta  (P)  -  Sideslip  angle 
CA  -  California 

Cdo  -  Coefficient  of  drag  at  a  =  0 

Centurion  -  Test  Pilot  School  student  in  Class  03A,  The  Wright  Stuff,  celebrating  100 
years  of  powered  flight. 

CG  -  Aircraft  Center  of  Gravity,  measured  in  %MAC 

Clo  -  Lift  curve  slope 

Cip  -  Rolling  moment  coefficient  due  to  P 

Ci6a  -  Rolling  moment  coefficient  due  to  aileron  deflection 

Ci6r  -  Rolling  moment  coefficient  due  to  rudder  deflection 

Clo  -  Coefficient  of  lift  at  a  =  0 

Cip  -  Rolling  moment  coefficient  due  to  rolling  moment 

Cmo  -  Pitching  moment  coefficient  due  to  a 

CmSe  -  Pitching  moment  coefficient  due  to  elevator  deflection 

Cmo  -  Pitching  moment  coefficient  at  a  =  0 

Cmq  -  Pitching  moment  coefficient  due  to  pitch  rate 

Cnp  -  Yawing  moment  coefficient  due  to  p 

Cn8a  -  Yawing  moment  coefficient  due  to  aileron  deflection 

Cn6r  -  Yawing  moment  coefficient  due  to  rudder  deflection 

Cnr  -  Yawing  moment  coefficient  due  to  yaw  rate 

CRP  -  Cramer-Rao  parameter 

Delta  (8)  -  Generic  flight  control  surface  deflection 

Delta  a  (6a)  -  Aileron  Control  Surface  Deflection 

Delta  e  (8e)  —  Elevator  Control  Surface  Deflection 

Delta  r  (8r)  -  Rudder  Control  Surface  Deflection 

DSN  -  Defense  Services  Network 

EOM  -  Equations  Of  Motion 

FTE  -  Flight  Test  Engineer 

FTN  -  Flight  Test  Navigator 

FTT  -  Flight  Test  Technique 

Gamma  (y)  -  Flight  path  angle,  climb/dive  angle 

h  -  height 

HAVE  TRIM  -  Project  title  bestowed  upon  this  testing  by  the  test  team. 

JON  -  Job  Order  Number 

JP-5/JP-8- Jet  Fuel 

LJ-25  -  Leaijet  25,  the  test  aircraft 

M  -  Mach  number 

MAC  -  Mean  Aerodynamic  Chord 

MOP(s)  -  Measure(s)  Of  Performance 


169 


DECEMBER  2003 


Edwards  Air  Force  Base 
Air  Force  Flight  Test  Center 


n*  -  Load  Factor  in  x-body  axis 
ny  -  Load  Factor  in  y-body  axis 
Oz  -  Load  Factor  in  z-body  axis 
P  -  Aircraft  Roll  Rate 
PRR  —  Preliminary  Report  of  Results 
OH  -  Ohio 

PA  -  Pressure  Altitude 

PEST  -  Parameter  Estimation  software 

Phi  ((]))  -  Euler  Bank  Angle 

PRR  -  Preliminary  Report  of  Results 

PTI  -  Programmed  Test  Input 

q  -  Dynamic  Pressure 

Q  -  Aircraft  Pitch  Rate 

R  -  Aircraft  Yaw  Rate 

R1  -  Region  /  Regime  1 

R2  -  Region  /  Regime  2 

RTO  -  Responsible  Test  Organization 

SHSS  -  Steady-Heading  Sideslip 

SLUT  -  Straight  and  Level,  Unaccelerated  Flight 

SPORT  -  Space  Positioning  and  Orientation  Radar  Tracking 

sqrt  -  Square  root 

T  -  Thrust 

Theta  (0)  -  Euler  Pitch  Angle 

TIM  -  Technical  Information  Memorandum 

TPS  -  Test  Pilot  School 

TPS/DOR  -  Test  Pilot  School  Operations  Desk 

UAV  -  Unmanned  Aerial  Vehicle 

V  -  Volts 

V  -  Velocity 

VDC  -  Volts,  Direct  Current 

VMC  -  Visual  Meteorological  Conditions 

VSS  -  Variable  Stability  System  installed  on  the  test  aircraft 

Vstaii  -  Stall  Airspeed 

W  -  Aircraft  Weight 

WUT  -  Wind-up  Turn 

lOCRP  -  Ten  times  Cramer-Rao  parameter 

95  SPTG/LGS  -  95*'’  Logistics  Squadron,  Edwards  AFB 
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